402  788  ~  *OT!'  40^788 


JL3.J 


1  MARCH  1963 


^epont: 

FOUNDATIONAL  RESEARCH  PROJECTS 


OCTOBER-DECEMBER  1962 


NAVAL  ORDNANCE  LABORATORY  CORONA 

CORONA,  CALIFORNIA 


NAVWEPS  REPORT  8141 


FOR  ERRATA 


AD _ 

402  788 


THE  FOLLOWING  PAGES  ARE  CHANGES 


TO  BASIC  DOCUMENT 


Z&IxoJt  <ry 


5  June  1963 


ERRATUM 


QUARTERLY  REPORT:  FOUNDATIONAL  RESEARCH  PROJECTS, 
OCTODER-DECEMBER  1962,  NAVWEPS  Report  8141,  U.  S.  Naval 
Ordnance  Laboratory,  Corona,  California,  1  March  1963. 

On  page  68,  immediately  preceding  the  last  paragraph,  the  line  should 
read: 

b  =  0.387  millivolts 


NAVAL  ORDNANCE  LABORATORY  CORONA 


W.  R.  KURTZ.  Ctfl.i  U8N 
Commanding  Offlctr 


F.  S.  ATCHISON.  Ph.  D. 
Tichalctl  Director 


FOREWORD 


This  report  includes  papers  on  various  aspects  of  work 
performed  during  the  second  quarter  of  fiscal  year  1963  under 
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Naval  Ordnance  Laboratory  Corona,  in  all  fields  of  science 
and  technology."  The  projects  are  supported  by  WepTask  R360 
FR- 104/21 1- l/ROi 1-01-001 .  Additional  support  is  provided 
the  High  Temperarure  Polymer  Program  by  WepTask  RMGA- 
41-031/21 1-1/F009-07-002. 


C.  J.  HUMPHREYS 
Head,  Research  Department 


CONTENTS 


Page 

FERROMAGNETIC  RESEARCH 

Spin  Wave  Damping  Constants  for  Disordered  Crystals, 


by  C.  A.  Roberts .  1 

Mean  Free  Path  of  Magnons,  by  C.  A.  Roberts . 20 

HIGH  TEMPERATURE  POLYMER  PROGRAM 

Synthetic  Studies,  by  D.  L.  Herring  and  C.  M.  Douglas . 21 


INFRARED  ATOMIC  SPECTRA 

Observation  and  Interpretation  of  4f-5g  Transitions  in  the 
Spectra  of  Alkalis  and  Noble  Gases,  by  C.  J.  Humphreys 


and  E.  Paul,  Jr . . 29 

NONAQUEOUS  ELECTROCHEMISTRY 

The  Electroreduction  of  Nitrobenzene  in  Liquid  Ammonia 

Solutions,  by  W.  S.  Harris . 47 

Ionic  Melt  Electrochemistry,  by  R,  E.  Panzer . 53 


NON  UNEAR  TRANSMISSION  ONES 

The  Semigeneralization  of  the  Transmission  Line  Equations, 

by  J.  R.  Aid  ay . 55 

SEMICONDUCTOR  MICROWAVE  MICROELECTRONICS 

Substrates,  Circuitry,  and  Semiconductors,  by  R.  W.  Yancey.  .  .  59 

SEMICONDUCTOR  PHYSICS 

Accurate  Measurements  of  a  Magnetic  Field  by  Means  of  a 

III-V  Hall  Effect  Device,  by  K.  I.  Crowder  and  H.  Filler  ...  67 

SOLID  STATE  SPECTROSCOPY 

The  Infrared  Absorption  Spectrum  of  Ce^+  in  LaClj,  by 


R.  A.  Buchanan,  H.  H.  Caspers,  and  J.  Murphy . 73 

Symmetry  Coordinates  for  LaClj,  by  J.  Murphy  and 

H.  H.  Caspers . 86 


SPECTRAL  EMITTANCE  OF  SOUDS 

Spectral  Emittance  Measurements  of  Some  Commercial 

Transparent  Solids,  by  D.  L.  Stierwalt,  D.  D.  Kirk,  and 
J.  B.  Bernstein . 


103 


FERROMAGNETIC  RESEARCH 


SPIN  WAVE  DAMPING  CONSTANTS 
FOR  DISORDERED  CRYSTALS 

by 

C.  A.  Roberts 


Because  the  damping  constants  for  spin  wave  excitation  are  mostly 
phenomenological  in  nature,  a  desirable  theoretical  objective  is  to  cal¬ 
culate  these  constants  in  terms  of  fundamental  constants  of  disordered 
crystals. 

The  excitation  of  spin  waves  in  ferromagnetic  films  has  been  amply 
demonstrated.*  However,  at  the  present  time  there  is  relatively  little 
understanding  of  the  splitting  of  the  resonance  peaks  and  the  nature  of 
the  damping.  The  purposes  of  the  calculations  presented  in  this  report 
are  to  determine  (1)  the  frequency  shift  of  the  normal  modes  (spin  wave 
modes),  and  (2)  the  damping  of  the  spin  wave  due  to  the  disorder  of  the 
sample. 

The  ideas  used  in  these  calculations  are  based  on  the  recent  work 
of  Maradudin  et  al.^  who  calculated  the  frequency  spectra  for  a  dis¬ 
ordered  lattice.  Their  model  is  a  one-dimensional  lattice  of  N  atoms, 
in  which  we  consider  N  -  n  of  the  atoms  to  be  of  mass  m  and  n  atoms 
of  mass  M.  We  now  consider  the  m-atoms  to  be  the  host  lattice  and 
the  M-atoms  to  be  distributed  at  random.  If  we  define  a  mean  mass  as 

Mt  =  Pm  +  (1  -  P)M 

then  we  can  consider  an  average  mass  M  perturbed  by  the  probability 
P  that  the  fth  atom  is  m  and  the  probability  1  -  P  that  the  f  +  k  atom 
is  M. 


1C.  Kittel,  Phys.  Rev.,  Vol.  110,  p.  1295(1958). 

2 

A.  A.  Maradudin,  G.  H.  Weiss,  and  D.  W.  Jepsen,  J.  Math.  Phys., 
Vol.  2,  pp.  349-369  (1960). 
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In  the  caee  of  spin  waves  in  a  crystal,  there  is  a  similar  disordered 
situation.  Here  the  characteristic  that  is  random  is  either  the  exchange 
constant  or  the  z  component  of  the  spin. 

Let  us  examine  the  microscopic  equations  of  motion  for  a  one¬ 
dimensional  spin  wave. 


where 


H.  =  2J  £g\g\ 

1  *7.  i  j 

i<J 

Using  the  semiclassical  approximations,^  we  have 

■fiST  =  2j[  5T  x  (6ST  +  a2V2?)]  (1) 

where  a  is  interatomic  spacing  and  where  the  quantity  in  parentheses  is 
Heff»  the  effective  field  at  the  j  atom  due  to  all  other  spins  in  the  struc¬ 
ture.  Since  Jf  x  6JT  =  0,  equation  (1)  reduces  to 

ST  =  [  ?  X  V2?]  (2) 


In  component  form,  we  have 


Sz 


is  assumed  to 
The  equation 


S 

x 


s  v2s 

z  y 


S 

y 


s  v2s 

z  x 


S  3  0 


z 

be  constant  in  space-time  and  equal  to  Sq. 
for  Sx  is 


C.  Kittel,  Solid  State  Physics,  New  York:  John  Wiley  and  Sons, 
Inc.  (1956). 


2„  v  2 


v4s 


(3) 


For  the  continuous  case  in  which  Sx  ~  e^(ut"^x),  the  equation  is 


or 


2g  /2I*2Sofk4g 

■“  So  *  -l— S— J  k  So 

X  '  '  X 


,  2Ja2S  , 

2  o,2 

CO  =  - - k 


(4) 


where  k  =  2n/X . 

The  dispersion  relation  for  long  wavelength  spin  waves  is  given  by 


E  =  -hu  =  2Ja2S  k2 
o 


The  finite  difference  equation  for  Sx  is 


*  «s] 


for 


S  ~  S  si<“‘-k<»l  =  s  (t)e-ik,a 

XX  X 

to  o 

i  t 

where  kQa  =  2irik/N.  For  periodic  boundary  conditions,  we  have 


and  as  the  general  solution 


S  =  S 

xN+i  x< 


k=*f 


1  \  '  ^  /2frikf  \ 

•  L  <Vx»(— ) 


k=-T+1 


(5) 


(6) 


(7) 


(8) 
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where  N  is  the  number  of  atoms  between  periodic  boundaries.  Thus 


Since 


and  Qk  =  Q*r 


Tj£exp[¥l(k  '  =  4(k  •  k'>  '  {J  k  5  k' 

l 

TjSy-S.  =  V  £££QtkQke’‘p[ir<(k  •  k'']=  EK 

s  a  I  «  t  i  <  . 


(9) 


k  k'  t 


Upon  substituting  equation  (8)  into  equation  (5),  the  finite  difference 
equation  for  the  spins,  we  get 


~p(^)=-(v)E^M^ 


+  exp 


4  exp 


4exp[~~  (f  +  1)J 


(10) 


For  each  k  we  then  have 

2 

/2JS  \ 

°k =  -2nr)  Qk(co*-^  -  4co,^h  4  3) 

and  with  trigonometric  substitutions 

X  .  So\  /,  2trk  ,  2  2irk\_ 

°k  +  4\  S  '/  (1  -  2co.— tco.  — JC^  =  0 

Note  that  for  Zwk/N  <  lor  kQa  <  1,  equation  (12)  reduce,  to  , 

and  the  normal  expression  for  long  wavelength  spin  waves  becomes 

..  I™  of  4 

+  l  ^  1  (kfc)  Qj^  =  0.  The  dispersion  relation  in  this  motion,  and  for 
subsequent  work,  is  then 


(11) 


(12) 


(13) 


where 


(14) 


Let  us  now  consider  the  effect  of  perturbing  the  normal  order  of  a  homo¬ 
geneous  line  of  spins.  By  this  is  meant  that  the  environment  of  a  par¬ 
ticular  lattice  site  is  changed  as  one  proceeds  down  the  line. 

The  average  exchange  is  defined  as 

J  =  PJj  +  (1  -  P)J2  (15) 


which  means  that  the  environment  of  atom  1  has  a  probability  P  and  that 
of  atom  2  has  a  probability  1  -  P.  The  equation  of  motion  for  the  system 
is  then 


'(i?)2<J' '  J,2[s'*2  * s'-2 '  4(s'+i  +  s'-*)+ 6S'J 


(16) 


where  Jj  is  the  exchange  constant  of  the  Ith  atom.  When  Jj  =  J,  the 
normal  equations  for  a  homogeneous  line  are  recovered.  Upon  substi¬ 
tuting  the  normal  modes  into  equation  (16),  we  have 


/2iTikf\  v— *  /2irik!\ 

Qk  exK— )  =  ■  X,  “k  Qk '‘i5  (— ) 


(17) 


After  multiplying  both  sides  of  equation  (17)  by  ^j^exp  (~ * )  an<l 

k' 

then  summing  over  I,  we  have 
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k*  i 


'xp[2”i<kNh'M]  -  -jjJ/k  Q‘‘“P 


ki 


By  using  the  orthogonality  property  of  the  exponentials  of  equation  (9). 
the  resulting  equation  for  the  Q^i  is 


(18) 


In  the  more  compact  form,  we  have 

^k  +  “k  Qk  =  "52®kk'Qk' 
k' 


(19) 


where 

•kk.  ■  -  *)  “*[¥<“  •  k,)<]  *  0  (20> 

1=2 

if  all  3f  =  J.  Here  4^)  can  be  considered  as  a  negative  pumping  term 
of  the  type  that  is  found  in  parametric  amplifier  equations.  In  the  spe¬ 
cial  case  in  which  all  =  0  except  and  =  *0cos2u>t,  equa¬ 

tion  (19)  takes  the  form 


•  •  2 

Q  +  (w  -  ®  cos  2ut)Q  =  0 
o 

This  equation  has  the  approximate  solution 


7  *0  \ 

[  sin  wt  -  cos  wt  + - r-Msin  3wt  -  cos  3wt) 

A  16  u  / 

Thus  Q  is  damped  by  #0. 


(21) 


(22) 


The  following  method  for  the  solution  of  equation  (19)  is  the  same 
as  that  used  by  Maradudin  et  al.,  ^  with  the  exception  of  some  small 


modifications.  Since  their  method  is  not  widely  used,  it  is  discussed  at 
length  below. 

The  transform  of  Qk(t)  is 


-St. 


c(s)  £  e  ““Q^dt 


(23) 


With  this  expression,  equation  (19)  takes  the  form 

SQk(°)  +  4(0)  ,  _ 

Vs,=  - 

k  kj 


where 


.  J-  2 

d,  =  S  +  u, 
k  k 


This,  of  course,  is  still  an  integral  equation  in  qk^(S). 
Upon  iterating  equation  (24),  we  have 

,k(S)  =  V(S)ck(S)  +  EDkn,'S»Cm'S» 

m^k 

where 


ck(S)  =  SQk(0)  +  Qk(0) 


and 


_  km  1  . _ 

DU„(S)  =  ; - — r— f® 


$  9 

kkj  k^m 


km  d,  d.  d  I  km  -  -  „ 

k  k  m  \  k, 

kl  1 


•‘k/k.kzV 

dk  dk 

kjk2  kl  k2 


(24) 


(25) 


(26) 
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Next,  the  authors  define  the  operator  Gj^S)  as 
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°k«s>  *  *kk  -  E 

ki 


a  9 

•  k'k,  k,k 


1  1  +  jV 


‘1 


klk2 


^kk^kjk^^k^k 

K1  *2 


(27) 


where  the  prime  means  that  no  summation  index  kj  shall  equal  k. 
The  relationship  of  to  G^  is  simply 

°k  =  3i?"5^ 

so  that  the  solution  to  the  inverse  transformation  is 


(28) 


Gk(S) 

VS)  =  rV + 


ZGm(s) 
m 

d  +  G 
x  y.  k  k 
m^k 


(29) 


where  only  the  diagonal  part  of  the  operator  Dj^m  is  used. 

The  ensemble  average  of  qfc(S),  the  average  over  all  atom  sites,  is 


and 


/CAS)  \  /G  (S)  \ 

<q(sl> = W  ^  (w 

‘  dk  *  <V 


(30) 


(31) 


The  solution  is  then  just  the  inverse  transformation  to  (Qfcft)). 

The  reason  for  using  only  diagonal  terms  is  that  the  average  is  per¬ 
formed  only  over  ^Qj<(t)>  and  not  over  ^Qj^Jt)  Qjj(t)^ .  In  the  latter  expres¬ 
sion,  the  diagonal  terms  do  contribute. 


Thus  we  find  that 


(32) 


or 
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(33) 


The  second  term  vanishes  for  m  /  k. 

The  contour  for  equation  (32)  is  the  usual  Bromwich  contour  for 
inverse  Laplace  transforms,  as  shown  below. 


The  behavior  of  ^Qk(t)^  is  governed  principally  by  the  denominators  of 
equation  (33).  The  expression 

N  *kk  *k  k 

°k<s)*»kk-  £  >  z  forS>0 

k,As  *“k 

is  the  first  approximation  to  equation  (27)  and  has  a  finite  number  of 
poles.  Since  +  d^  +  will  have  a  finite  number  of  zeros,  the 

perturbed  system  will  tend  to  oscillate.  For  an  infinite  number  of  par¬ 
ticles,  the  pole  density  becomes  infinite  in  a  finite  length  (um^  -  wm). 

If  we  substitute  iy  for  S, 


Gk^>  =  *kk  - 


There  are  2  N  poles  between  w  -  w  at  y  =  ±w  ±  w. 


N 


$  $ 
kkj  kjk 


(34) 
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Wn 


Ai  can  be  seen  from  the  plot  of  equation  (34),  there  is  an  infinite  num¬ 
ber  of  poles  a»  N  ■*  «. 


v  -  v =  +  co,2'~<V 11 


(35) 


-  COS 


2ir(k  +  1)  ,  .  2  2 irk  2irkl 

~  n  • 1  +  c°’  —  +  co,TrJ 


From  equation  (35)  we  obtain 


*  fv  T  f  2  2ir ..  .  2  2ir  2w(k  +1)  2ir .  1 

=  V 4yJ  jcos  —  (k  +  1)  -  cos  — k  +  cos — 1  -  cos— kj 

For  N  large,  this  becomes 

■ 


Aw, 


for  2k  +  1 


and,  as  N  -*  °°  and  2irkja  —  0 


Aw, 


*  ) 


for  fixed  kna  <  1 


(36) 


The  expression  for  ^Qfc(t)^  is 


St 


QJc)  /,cti“  -  s 

<QJ')>  ■  -sr  /  i — r - 

Jc-ioo  S  +  wk  +  <Gk(S)> 


dS 


Q ,(0)  fcHao  St 

k  /  e  c 

2iri  J  .  2  2 

•/c-i«  S  +  w,  + 


(37) 


dS _ 

<GJS)> 
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Now  consider  the  integral 


^.c+ioo 

&[  - 

/  2  2 
•'c-i*>  S  +  u,  + 


eStdS 


9  9 

kk.  kk 


.  4 

k  kk 


^  1 


which,  on  passing  to  the  limit,  becomes 


_1_  /•ct‘“ _ e^dS _ 

2ffi  /  2  2  if 

J  .  S  +  w  (k)  +  9..  -  ■= —  /  — 5 - dk 

-Vi”  kk  2V  s2 *  atkj)  1 


Aa  S  -*y  and  c  -*■  0,  we  have 


■  /*c+i«o  St 

1  /  _ e  dS _ 

2-rri  /  ?  ,  .  /K(kk  )dk 

y.-«-  -y* -*«♦•« -hj-T-T-\ 

J  w  (k  j)  -  y 


(38) 


(39) 


(40) 


This  expression  reveals  the  fact  that  the  integral  does  not  have  a  simple 
pole  and  so  the  integral  in  the  denominator  prevents  the  integrand  from 
being  analytic  across  the  imaginary  axis. 

The  cut  is  introduced  in  the  S-plane  from  u,,,  to  -wm. 
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Thu* 


where 


lim  G  (iy  ±  t,)  =  -  K  (y)  *  U.  (y) 

t!  — 0+  K  K  K 


-  ■  *kk  + 


ft  ft 

,  kk.k.k 
^  1  1 

;R-y2), 


(41) 


(42) 


and 


Vyl  *  ’"ir,yS*kk1®k1k6(“k1  - y2) 


(43) 


Upon  changing  variables,  s  =  iy  +  c  as  c  -*•  0+,  and  we  obtain 


V4* 


i  r°° 

'  2ir  /  2 

-/.«  -y  +  u. 


k  <Ky(v)>  +  i<vk(y)> 


-dy 


(44) 


The  denominator  of  the  integral  is  juBt  the  average  of  the  operator 

D^(s)  =  ~2  2 


s  +  u  +  G, 


which  can  be  written  as 


Dk(»)  =  [t  -  z  -  g(z)] 

and  which  has  dissipative  properties. * 
Equation  (44)  can  be  written  as 

+°°  iyt 


V*> a  zi 


l  r  _ e_] 

ukly)  + 


iVk*y) 


dy 


(45) 


N.  M.  Hugenholtz,  The  Many  Body  Problem,  New  York:  John 
Wiley  and  Sons,  Inc.,  pp,  33-35  (1959). 


where 

<vk(y)>  =  0  for  |y|  >  yQ  Y0  =  |wn|  (46) 

that  is,  the  imaginary  part  of  the  denominator  vanishes  for  regions 
beyond  the  maximum  frequency.  Also, 

<uk(y)>  =  -y2  +  wk2  +  =  0  for  y  >  yQ  (47) 


If  the  probability  value  P  for  a  lattice  site  is  equal  to  1  or  0,  then 
Vk(y)  =  0  for  all  y  and  Uk(y)  reduces  to  -y2  +  u>j^  =  0,  and  the  normal 
unperturbed  frequencies  are  recovered. 


It  is  the  finite  and  nonvanishing  part  of  ^Vk(y)^  that  gives  rise  to 
dissipation.  This  is  dissipation  in  the  following  sense:  The  state  of  the 
system  changes  in  time  r  irreversibly  goes  into  a  more  complicated 
configuration,  but  there  is  no  violation  of  the  conservation  of  energy. 

The  argument  of  Hugenholtz4  can  be  used  to  estimate  the  long  time 
behavior  of  the  integral.  It  must  be  emphasized  that  since  his  method 
would  give  a  vanishing  result  for  the  normal  modes,  it  would  not  give 
the  complete  time  dependence  of  the  system. 


If  we  expand  Uk(y)  +  iVk(y)  about  the  points  y  =  +yk  and  y  =  -yk  and 
note  that  Uk(yk)  =  0, 


■*L  -’'te 


1 


(yj  +  (y  -  yt/U.((y. )  +  iv  (y) 


k'  k'  k 
1 


(48) 


V'V +  (y  +  yk,uk'('yk' +  ivk('yk’ 


] 


dy 


for  4k  >  0.  Combining  these  terms  with  the  given  notations 


we  have 


and  Fk  =  NkVk^k> 


uk(y)  =  %2) 


and  Uk,(-y)  =  “^.(y) 


which  may  be  plotted  as 


(49) 


(50) 


13 


Thus 


V*> 


_ 1_ 

y  +  yk- 


dy 


(51) 


Ik(t)  »  2Nke  k  Binykt  (52) 

The  perturbed  spin  wave  mode  is  then 

r  "rkl  i 

<Qk(t)>  ~  2Nk  [Qk(0)e  ykcosykt  +  Qk(0)e  sinyktj  (53) 

Upon  substitution  of  the  above  equation  in  the  expression  (8)  for  the  nor¬ 
mal  modes,  we  have 


k=  +N/2 


-r,  t 


VTn  E  2Nk6k(0).  .i»yk,.xp(^)  ,54) 

1  k=-(N/2)+l 

If  we  assume  Qk(0)  ^  0  and  Qk(0)  =  0,  then 

%  ■  toE2nA<°>  «■>  [*(v  -  ^)]  <55> 


This  is  the  general  solution  to  spin  waves  propagating  in  a  disordered 
medium,  where 

2NkQfe(0)  =  Ak  amplitude 

rk  =  damping  constant 

and  yk  =  the  shifted  normal  mode  frequency,  which  is 

determined  from 
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(56) 


+  a>k2  +  <Gk(iy  +  »1)>  as 


Tl 


Physically,  the  disordered  spin  wave  medium  damps  a  normally  propa¬ 
gating  spin  wave  for  long  periods  of  time. 


It  must  be  emphasised  that  this  is  not  the  only  mechanism  that 
brings  about  a  finite  amplitude  to  a  spin  wave  resonance. 


The  ensemble  average  of 


jA(k  -  kj) 


The  average  for  (&kki®kjk?^  is  calculated  in  the  same  manner  as  that 
used  by  Maradudin  et  al.f  tnat  is,  we  separate  the  diagonal  terms  from 
the  nondiagonal  terms  and  sum  and  average  separately. 


kiV’  = 


w  k  w  k? 


•  exp 


2ttL 

N 


[<* 


+ V 


2,  K 

w  kiw  k2 


N 


(Contd. ) 
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(59) 


S4HB-0 


exp  {"^[(k  *  kl)j  +  (kl  ‘  k2)i]} 


We  now  average  each  term  separately. 

2.  2, 

o 

<* 


,*klk2> — (((r  - ')  )Z »p[-¥(k  •  k2»] 

^  (60 

*  {Jt  ‘  ‘ )  ^  £  «P  {-^[<k  -  kl»  +  <kl  -  k2>']}^ 


Next,  we  add  and  subtract  the  diagonal  term 

2  2  2 
“*1  W  k2  fP2  ’  P1 


<VVk> *  ^i(k  -  ki»]  •  a,k  •  k2»  <61> 


‘1  12 


where 


+  (1 


+  0  - 


With  these  results,  we  can  now  find  explicit  expressions  for  the  damping 
constants. 

The  solution  for  the'  shift  of  the  normal  mode  frequency  is  reached 
in  the  following  manner:  From  equation  (56),  with  =  0,  we  have 


X*kkJ*k|k)> 


^  (“  > «,  -  y\ 


where  the  averages  may  be  obtained  from  equations  (58)  and  (61).  Then 


2 


Treating  the  sum  as  small,  we  have  to  the  first  order 


2 


We  may  proceed  in  the  manner  of  Maradudin  et  al.  and  calculate  the 
shift  in  the  aero  point  energy,  but  that  is  not  of  interest  here. 


The  damping  constant  was  given  in  equation  (49)  as 

rk  ■  VA> 

where  y^  is  the  shifted  normal  mode  frequency  and,  from  equation  (43), 

W  =<W>  ■  <’8in'r2Xk*k1k6K2  • 

k,  *  > 

Using  the  specific  results  for  the  spin  wave  problem  as  given  by  equa¬ 
tion  (61)  for  we  have 

W  ’  •"k2(l‘2  ‘  1‘i2)ykSE“k16(“k1  ’  yk2) 

=  -  co,ir)6("Ji  -  yk2) 

Using  the  relation 
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we  obtain 


where  Q  -*  irk/N.  We  let 

2  2  2 
X  =  u  (1  -  cos  2Q) 
m 


and  solve  for  cos  2Q.  If  we  use  the  fact  that  0  <  £  u»m  and  the  relation 


then 


6<X2  -  y2)  =  ^[6(X  -  y)  +  6(X  +  y)] 


=  0  y  >  u 

m 

The  other  function  needed  for  the  damping  constant  is 

*  |  V<V 

Since 

we  can  approximate  I  Uj^ ,  (y)  | : 
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|uk«  <y>|  “  | -  2yk  ■  <Ky. (y)>|  =-  2wk 

to  the  first  order  for  yk  cm  o»k.  By  combining  the  terms,  the  damping 
constant  for  spin  waves  as  a  function  of  wave  number  is 

*tV(l‘2  '  **1Z)  “k 

4ir  u 

m 


r  =t(p2  -  p,2)« 

m  4  v  2  1  /  m 


This  is  to  be  compared  with  the  case  of  phonons  where 


r  =  * 

m 


For  long  wavelength  spin  waves,  with  ka  «  1 


rk='?(,‘2-  l*l2)“m"  '  co,ka) 


or 


and  for  =  wm 


r 


k 


2 


u 

m 


2JS 

o 

h 


~  2  x 


-1 

sec 


^2  "  *1 


P(1 .  P,[(i* .  .)2  ♦  (£ 


where  Jj  and  J2  are  the  exchange  constants  for  each  component.  Thus 


rk=io,3P(i-P)[(ii-  +  ')2]<ka>3 

for  long  wavelength  spin  waves. 

Note  that  spin  wave  energy  is  proportional  to  (ka)2  and  spin  wave 
velocity  is  proportional  to  ka,  therefore  the  rate  of  energy  transfer  is 


1« 


1 


proportional  to  (ka)3.  For  example,  if  X  ~  100  A,  a  ~  1  A  and 
rk  ~  109  see" 


MEAN  FREE  PATH  OF  MAGNONS 
by 

C.  A.  Roberts 

The  lifetime  of  the  magnon  is 

4u >  «/ 1  -  (1  -  w  /u  )* 
l  m  km 


Tk  =  rk 


(^1  -  ^I2) 


The  mean  free  path  as  a  function  of  number  is 


du 


I  =  V  t  where  V  =  —  =  group  velocity 

K  K  £  QK 


Thus 


=  ,inka^‘  -  <*  - 


Mi  -  M2 

for  spin  waves  with  ka  «  1,  and  for  ka  =  1/100 


i,  =  1.6  x  10  cm 
k 
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HIGH  TEMPERATURE  POLYMER  PROGRAM 


SYNTHETIC  STUDIES 
by 

D.  L.  Herring  and  C.  M.  Douglas 
It  was  initially  postulated  by  the  authors  that 

r  NHz 

h 

L(C6H5)2P-N 

(1) 

would  react  with  the  three  classes  of  pentavalent  phosphorus  halides 
according  to  the  following  equations: 


Compound  (I)  +  R^PC^ 

- ►  HC1  +  (R  PN) 

L  J  |  X 

(1) 

Compound  (I)  +  PCI  - 

b 

-  hci  +  (R4ci2p3n3)x 

(2) 

Compound  (I)  +  RPCl^  — 

— ►  HCI  +  (R5C1P3N3)x 

(3) 

NH_ 


-----  P(C 


>H5>jC1 


and  yield  the  corresponding  cyclic  trimer,  along  with  higher  polymers. 
The  preparation  of  [  (CfcHsJ^PN^  by  the  equation  (1)  method  has  been 
described.  1  This  report  gives  an  account  of  the  preparation  of 
[(CfcHs^Pj^C^]  by  means  of  equation  (2). 


When  compound  (I)  reacted  with  PCI5  under  vacuum  over  the  range 
135-145*C,  about  75%  of  the  amount  of  HC1  was  evolved  that  would  be 
expected  from  the  reaction 


NH. 


NH. 


2  |  2 

L<C6H5>2P-N=P<C6H5>d 


Cl  +  PC15  - ►  [(C6H5)4P3N3C12]  +  4  HC1 


NAVWEPS  Report  7229.  Quarterly  Report;  Foundational  Research 
Projects,  April-June  1962,  pp.  19-31  (August  1962). 


From  the  product  mixture,  tome  (PNCl2)3  was  isolated,  together  with 
a  24%  yield  of  a  white  crystalline  compound  (m.p.  142-143*C,  lit. 
m.p.  137-140*0),  which  was  identified  on  the  basis  of  its  infrared 
spectrum  (Figure  1)  and  an  elemental  analysis  as  the  tetraphenyl  gemi- 
nal  dichlorophosphonitrilic  cyclic  trimer^ 


«5C6V  /C6H5 

/P\ 

N  N 

H5C64  |/C1 
H5C6A  /  "C1 

N 


(id 


The  isolation  of  the  geminal  dichloride  caused  considerable  diffi¬ 
culty  because  of  its  extreme  solubility  in  all  the  solvents  tried.  How¬ 
ever,  this  route  to  the  geminal  dichloride  appears  to  be  an  improve¬ 
ment  over  the  Friedel-Crafts  method,  ^  which  required  a  reaction  period 
of  six  weeks  and  gave  lower  yields. 

In  a  second  experiment,  compound  (I)  reacted  with  PCI5  in  CH3CN 
and  again  yielded  the  geminal  dichloride  (II),  but  the  yield  was  less  than 
the  amount  isolated  in  the  first  experiment  described  above. 

In  a  preliminary  experiment,  compound  (I)  reacted  with  C5H5PCI4 
under  vacuum  at  140-145*C  and  gave  a  new  [(CfcHs^PN]  x  compound 
(m.p.  163-164*C),  together  with  a  polymeric  material  which  analyzed 
only  approximately  for  [(CfcH5)5P3N3Cl) . 

In  a  series  of  four  experiments  in  which  reaction  times  were  varied 
from  45  to  90  hours,  (CfcHg^PCl  reacted  with  NH2NH2'HC1  in  refluxing 
tetrachloroethane^  and  yielded  [(CfcHs^PN^  and  a  compound  that  was 
tentatively  identified  as 


NH  Cl 

»  I 

<C6H5>2P-N  =  P<C6H5>2 

(III) 


2 

R.  A.  Shaw  and  F.  B.  G.  Wells,  Chem.  and  Ind.,  pp.  1189-1190 
(1959). 

3 

H.  Sisler,  H.  Ahuja,  and  N.  Smith,  Inorganic  Chemistry,  Vol.  1, 
pp.  84-88  (1962). 
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FIGURE  1.  Infrared  Spectrum  of  the  Tetraphenyl  Geminal  Dichlorophosphonitrilic  Cyclic  Trimer 


Pyrolysis  of  compound  (III)  under  vacuum  over  the  range  240-250*C, 
which  may  be  described  according  to  the  equation: 


NH  Cl 

II  I 

3(H5C6)2P-N  =  P(C6H5)2 


♦  3  HC1  +  2  [(C6H6)2PN]  3 


yielded  [  (CfcHs^PN]  3. 


EXPERIMENTAL 

Reaction  of  Compound  (I)  with  PCI5  Under  Vacuum 


In  the  drybox,  5.71  g  (0.0126  mole)  of  compound  (I)  was  mixed  with 
2.63  g  (0.0126  mole)  of  PCI5.  and  the  mixture  was  introduced  into  a 
heavy-wall  glass  ampoule.  The  ampoule  was  removed  from  the  drybox, 
attached  to  the  vacuum  line,  and  gradually  heated  in  an  oil  bath.  The 
reaction  mixture  melted  at  bath  temperature,  110-120*C,  and  evolution 
of  hydrogen  chloride  commenced  at  135-140’C.  The  reaction  was  con¬ 
tinued  for  8  hours  at  140-145*C,  with  continuous  removal  of  the  HC1. 
The  yield  of  HC1  was  830  ml/STP,  or  74%  of  the  amount  anticipated  for 
the  reaction 


r  NH,  NH,  -1 

II  2  I  2 

_<C6H5,2P-N  =  P(C6H5)2jC1  + 


PC15 - ►[(C6H5)4P3N3C12]  +  4HC1 


The  ampoule  was  detached  from  the  vacuum  line  and  opened  to  the  air. 
A  small  amount  of  white  crystalline  material  which  had  sublimed  into 
the  neck  of  the  ampoule  was  identified  by  its  melting  point  (109-110*C, 
Fisher-Johns  m.p.  block)  and  infrared  spectrum  as  (PNCl2)3- 


The  solid  material  remaining  in  the  ampoule  was  extracted  with  boil¬ 
ing  benzene.  A  portion  of  the  reaction  mixture  was  benzene-insoluble 
and  was  removed  from  the  solution  by  filtration.  The  insoluble  fraction 
was  extracted  with  400  ml  of  hot  ligroine  (boiling  range  65-75*C).  When 
the  ligroine  solution  was  allowed  to  stand  undisturbed  for  several  days, 
crystals  were  formed  which  melted  sharply  at  171  *C.  No  other  analyt¬ 
ical  data  are  available  as  yet. 


The  benzene  was  removed  from  the  soluble  portion  of  the  reaction 
product.  The  light-tan  residue  was  extracted  repeatedly  with  ligroine 
(boiling  range  90-120’C).  Prolonged  fractional  crystallization  from 
the  ligroine  solution  yielded  1.51  g  (24%)  of  compound  (II). 
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Calcd.  for  C,  .H,nP,N-Cl,:  C,  56.1;  H,  3.9;  P,  18.1;  N.  8.2;  Cl.  13.8; 

24  20  3  3  2  _  i  *7  i^a*/** 

m.p.  137-140  C 

Found:  C.  56.2;  H,  4.0;  P,  18.2;  N.  8.3;  Cl.  13.8 
m.p.  142-143*C 

The  material  that  was  not  soluble  in  ligroine  softened  over  the  range 
80- 100*C. 

Found:  C,  54.8;  H,  4.9;  P.  18.3;  N,  8.4;  Cl,  10.4;  O.  2.3 
Reaction  of  Compound  (I)  with  PCI5  in  CH3CN 


In  the  drybox,  4.2  g  (0.009  mole)  of  compound  (I),  dissolved  in 
100  ml  of  acetonitrile,  reacted  with  1.9  g  (0.009  mole)  of  freshly  sub¬ 
limed  PCI5  in  the  manner  previously  described. ^  One  of  the  products 
isolated  in  this  experiment  was  a  material  which  was  recrystallized 
from  ligroine  (boiling  range  90-120*C).  Its  melting  points  were  137- 
139*C  (Fisher- Johns  m.p.  block)  and  142-143*C  (sealed  tube).  This 
material  was  identified  as  the  geminally  substituted  dichlorotetraphenyl 
phosphonitrilic  trimer. 

Calcd.  for  C  H  p  N  Cl,:  C,  56.1;  H,  3.9;  P,  18.1;  N.  8.2;  Cl,  13.8 
24  20  3  3  2 

Found:  C,  56.3;  H,  3.9;  P,  18.1;  N,  8.4;  Cl,  13.6 
Preparation  of  Phenyltetrachlorophosphorane 

In  the  drybox,  50.4  g  (0.281  mole)  of  phenyldichlorophosphine  was 
dissolved  in  500  ml  CCI4  which  had  been  dried  by  distillation  from  P2O5. 
The  solution  was  introduced  into  a  1 -liter  stainless  steel  bomb  fitted 
with  a  valve  and  vacuum  line  connection.  The  bomb  was  attached  to  the 
vacuum  line,  cooled  to  -196*C,  and  charged  with  20  g  (0.281  mole)  of 
chlorine.  The  bomb  was  then  allowed  to  warm  to  room  temperature. 
After  72  hours  it  was  emptied  into  a  liter  flask  in  the  drybox.  Because 
of  the  great  solubility  of  C6H5PCI4  in  CCI4,  it  was  necessary  to  evapo¬ 
rate  the  solution  to  near  dryness  before  the  product  could  be  isolated 
(m.p.,  sealed  tube,  75-76'C;  lit.  75"C). 

Reaction  of  Compound  (I)  with  Phenyltetrachlorophosphorane 

In  the  drybox,  5.11  g  (0.011  mole)  of  compound  (I)  was  mixed  with 
2.82  g  (0.011  mole)  of  C6H5PCI4  and  introduced  into  a  heavy-wall  glass 
ampoule.  The  reaction  tube  was  attached  to  the  vacuum  line  and  gradu¬ 
ally  heated  to  175*C.  The  mixture  fused  at  90-95'C  and  began  to  evolve 


^NAVWEPS  7237,  Quarterly  Report:  Foundational  Research  Proj¬ 
ects,  July-September  1962,  p.  48  (December  1962). 
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HC1  at  120*C.  During  the  8>hour  heating  cycle  the  HC1  was  continuously 
removed  (total  amount  recovered,  720  mm/STP).  The  reaction  product 
Was  then  extracted  with  benzene  and  0.59  g  of  [(C{>H5)2PN]4  was  re¬ 
covered  by  fractional  crystallization.  The  benzene  was  evaporated  and 
the  crude  product  extracted  with  ligroine.  The  first  fraction  recovered 
from  the  ligroine  solution  was  a  small  amount  of  crystalline  material 
that  melted  sharply  at  163-164*C. 

Calcd.  for  f(C,H,),PN]  :  C.  72.3;  H,  5.1;  P,  15.6;  N,  7.0 
Found:  C,  72.1;  H,  5.1;  P,  15.6;  N,  6.8 

A  material  which  appeared  to  be  amorphous  (softening  point  95-100*C) 
Was  also  precipitated  from  the  ligroine  solution. 

Calcd.  for  (C.H^P^Cl:  C,  65.0;  H,  4.5;  P,  16.8;  N,  7.5;  Cl,  6.3 

Found:  C,  66.78;  H,  4.69;  P,  16.8;  N,  6.37;  Cl,  4.62 

Reaction  of  Diphenylchlorophosphine  with  Hydrazine 
Monohydrochloride  in  Tetrachloroethane 


Four  experiments  were  conducted  to  test  the  reaction  of  (C^Hs^PCl 
and  NH2NH2'HC1  in  tetrachloroethane.  Reaction  times  varied  from  45 
to  90  hours,  but  the  results  were  essentially  the  same.  The  following 
is  a  typical  experiment. 

To  21.0  g  (0.3  mole)  of  NH2NH2'HC1  slurried  in  200  ml  of  boiling 
tetrachloroethane  was  added  dropwise,  with  stirring,  a  solution  con¬ 
taining  66.3  g  (0.3  mole)  of  (C&H5)2PC1  in  200  ml  of  tetrachloroethane. 
After  the  mixture  was  stirred  at  reflux  temperature  for  72  hours,  a 
precipitate  of  NH2NH2'HC1  and  NH4CI  was  removed  by  filtration.  Next, 
the  excess  tetrachloroethane  was  removed  under  vacuum  and  an  amber- 
colored  mixture  of  solids  remained.  The  mixture  of  residual  solids 
was  extracted  with  1  liter  of  dry  acetone,  and  the  insoluble  material 
was  removed  by  filtration.  From  the  acetone  solution,  7.2  g  (11%)  of 
a  white  granular  mixture  (m.p.  217-230*C)  was  isolated  by  recrystal¬ 
lization.  The  material  was  probably  the  .desired  compound, 

NH  Cl 

II  I 

(CfcHs^P — N  =  P(C^H5)2>  contaminated  with  some  oxygen- containing 
impurity. 

Calcd.  for  C^H^P^Cl:  C,  66.6;  H,  4.9;  N,  6.4;  P,  14.2;  Cl,  7.9 

Found:  C,  68.8;  H,  5.1;  N,  6.2;  P,  13.2;  Cl,  5.1;  O,  1.6 


The  acetone-insoluble  material  was  then  extracted  with  1  liter  of 
boiling  benzene  and  the  insoluble  material  was  removed  by  filtration. 
From  the  resulting  solution,  by  the  process  of  fractional  recrystalli¬ 
zation,  3.2  g  (5%)  of  a  white  crystalline  material  (m.p.  267-269*C)  was 
isolated.  Concentration  of  the  remaining  benzene  solution  yielded 
10.3  g  (18%)  of  [(C6H5)2PN]3  (m.p.  224-226*C). 

When  a  sample  of  the  unknown  material  was  pyrolyzed  under  vacuum 
at  240-250*C,  95%  of  the  HC1  required  for  the  reaction 

NH  Cl 

II  I 

3(C6H5)2P—  N  =  P(C6H5)2 - ►2[(C6H5)2PN]3  +  3  HC1 

was  collected.  The  resulting  solid  product  was  identified  as  [(Cf>H5)2PN]3 
(m.p.  217-222*C).  The  conversion  of  this  material  into  [  (CfcHs^PN]  3 
in  the  above  pyrolysis  experiment  strongly  suggests  that  it  is  the  com- 
NH  Cl 

II  I 

pound  (C(,H5)2P — N  =  P(CfcH5)2-  However,  a  quantitative  analysis  of 
the  elements  was  only  approximate  for  this  compound. 

Calcd.  for  C26H2IP2N2C1:  C,  66.5;  H,  4.9;  P.  14.3;  N,  6.5;  Cl,  7.9 
Found:  C,  66.5;  H,  4.9;  P,  14.4;  N,  8.9;  Cl.  5.6 
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INFRARED  ATOMIC  SPECTRA 


OBSERVATION  AND  INTERPRETATION  OF  4f-5g  TRANSITIONS 
IN  THE  SPECTRA  OF  ALKALIS  AND  NOBLE  GASES 

by 

C.  J.  Humphreys  and  E.  Paul,  Jr. 


The  experimental  program  during  recent  months  haB  emphasized 
the  observation  of  selected  spectra  in  the  4-micron  region.  Owing  to 
repeated  detector  failures,  the  amount  of  data  accumulated  has  been 
relatively  small.  It  has  been  sufficient,  however,  to  demonstrate  the 
feasibility  of  observation  in  this  region.  Furthermore,  the  interpre¬ 
tation  of  new  data  on  the  spectra  of  the  noble  gases  has  demonstrated 
the  validity  of  the  treatment  of  one  aspect  of  the  theory  of  complex 
spectra,  something  heretofore  impossible  owing  to  lack  of  observations 
in  a  region  previously  unexplored  by  high-resolution  methods. 

Discussion  of  the  experiments  is  limited  to  a  few  comments,  since 
in  essentials  the  methods  remain  the  same  as  previously  reported. 
Microwave -excited  electrodeless  source  tubes  with  inblown  windows 
oriented  end -on  have  proved  satisfactory  for  the  region  covered.  The 
response  of  coolable  lead  sulfide  cells  mounted  in  miniature  Dewars  of 
a  variety  of  designs  is  found  not  to  fall  off  to  a  significant  extent  for 
wavelengths  short  of  4  microns.  The  location  of  the  practical  limit 
near  this  point  still  permits  observation  of  essentially  all  the  interest¬ 
ing  features  of  the  spectra  under  study.  A  much  more  serious  handi¬ 
cap  associated  with  the  use  of  these  detectors  is  a  tendency  of  the  De- 
wars  to  lose  vacuum  or  even  to  fracture,  something  apparently  associ¬ 
ated  with  the  cooling  and  heating  cycle.  Some  consideration  is  being 
given  to  developing  new  designs  for  the  Dewars  since  the  loss  of  time 
from  detector  failure  is  becoming  intolerable. 

It  has  also  been  necessary  to  replace  the  15,000  lines-per-inch 
grating  with  dne  of  7,500  lines  per  inch,  since  the  former  can  only 
reach  3.3  microns  in  the  first  order  in  the  Littrow  arrangement.  The 
relatively  small  dispersion  associated  with  practical  scanning  speeds 
is  not  so  much  of  a  handicap  as  it  might  appear  when  working  near  4 
microns.  With  reasonable  effort  the  chart  features  can  be  measured 
to  not  much  better  than  1  A.  It  is  noted,  however,  that  at  4  microns 
a  change  of  1  A  is  equivalent  to  a  wave-number  variation  of  0.06  cm"*. 
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In  an  earlier  period  of  radiometric  observation  of  line  emission 
spectra  connected  with  the  name  of  Paschen,  *  a  feature  was  found  in 
several  spectra,  notably  those  of  the  alkalis,  near  4  microns.  It  was 
interpreted  as  the  transition  between  the  unresolved  doublet  levels 
associated  with  the  f-  and  g-configurationa.  The  line,  shifting  slightly 
to  shorter  wavelengths  with  elements  of  increasing  atomic  number, 
was  very  near  to  that  associated  with  the  transition  between  the  fifth 
and  fourth  orbits  of  hydrogen,  which  was  discovered  by  Brackett2  more 
than  ten  years  later.  Illustrated  here  is  a  kind  of  degeneracy  related 
to  all  atomic  spectra,  resulting  from  the  essentially  hydrogenic  char¬ 
acter  of  levels  from  peripheral  configurations.  It  is  expected,  there¬ 
fore,  that  a  fairly  conspicuous  feature  may  be  observed  near  4  microns, 
of  a  complexity  dependent  on  the  electron  configuration.  The  precision 
of  Paschen's  observations  was  surprisingly  good  when  one  considers 
the  energy  limitation  of  his  thermal  detectors  and  the  fact  that  the  lines 
are  of  only  moderate  intensity. 

A  short  section  of  cesium  record,  shown  in  Figure  1,  illustrates 
the  feature  in  question,  another  Ritz  combination  in  Cs  I,  and  two  argon 
lines.  The  wavelength  of  the  Cs  I  line  of  leBser  intensity  is  calculable 
from  Kleiman'a-*  energy-level  values.  The  first-order  Ar  I  lines  are 
included  among  those  described  in  the  following  section  of  this  report. 

The  features  appearing  in  the  record  are; 


A 

<t  ,  cm'* 

Transition 

Cs  I 

39177.66 

2551.779 

5p67p2Po|  "  5p68s2S0| 

Cs  I 

39418 

2536.2 

5p64f2F«  .  5p65g2G 

Ar  I 

39659.75 

2520.76 

«[  1*1  1,2-  5«I  2*15.3 

Ar  I 

39899.65 

2505.60 

4f[«*U,5  -  5*l  5*  15.6 

Of  much  greater  theoretical  interest  are  the  features  associated  with 
the  transitions,  occurring  in  the  4-micron  region,  between  the  levels 
originating  respectively  in  the  configurations  p^f  and  p^g  of  the  noble 
atmospheric  gases.  These  configurations  are  analogous  with  pf  and  pg, 
which  are  characteristic  of  fourth-group  elements,  notably  germanium. 


F.  Paschen  and  R.  Gotze,  Seriengesetze  der  Linienspektren, 
Berlin:  Springer  (1922).  * 

2 

F.  S.  Brackett,  Astrophys.  J.,  Vol.  56,  p.  154  (1922). 

3 

H.  Kleiman,  J.  Opt.  Soc.  Am.,  Vol.  52,  p.  441  (1962). 
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which  has  been  the  subject  of  an  investigation  by  the  authors'*  in  coop¬ 
eration  with  K.  L.  Andrew  of  Purdue  University.  Application  of  the 
Pauli  exclusion  principle  to  configurations  involving  a  given  number, 
designated  n,  of  p-electrons  leads  to  the  same  number  of  levels  with 
the  same  set  of  J-values  as  that  obtained  for  (6  -  n)  p-electrons.  In 
general,  the  order  of  the  J-values  is  reversed  for  the  complementary 
cases.  When  the  outer  electron  assumes  large  f-values,  as  in  the 
instances  of  f  or  g,  the  distribution  and  quantum  number  characteriza¬ 
tion  of  the  levels  are  described  closely  by  the  pair-coupling  approxi¬ 
mation,  probably  first  formulated  by  Shortley  and  Fried-*  in  1938  for 
the  d^f  and  d^g  configurations  in  Cu  II  and  for  the  p^d  and  p^f  configu¬ 
rations  in  the  rare  gases.  A  general  treatment  of  intermediate  or  JJ 
coupling  was  given  by  Racah^  in  1942.  Our  own  recent  work  on  the 
experimental  determination  of  the  p^f  levels  of  Ar  1  is  described  in  an 
earlier  report^  in  which  the  physical  conditions  appropriate  to  the  appli¬ 
cation  of  the  Racah  theory  for  intermediate  or  J I  coupling  —  namely, 
electrostatic  interaction  relatively  weak  compared  with  the  spin-orbit 
interaction  of  the  parent  ion  but  strong  compared  with  the  spin  coupling 
of  the  external  electron — lead  to  two  groups  of  close  pairs  of  levels, 
the  groups  being  separated  by  the  interval  between  the  two  levels  of  the 
ion  limit.  In  the  extreme  case  where  the  spin  coupling  of  the  external 
electron  is  vanishingly  small,  the  components  of  the  close  pairs  coalesce 
so  that  instead  of  twelve  levels  for  the  configurations  under  discussion 
the  number  would  be  expected  to  be  reduced  to  six,  each  with  two  J- 
values.  Owing  to  a  degeneracy  that  appears  in  the  evaluation  of  the 
energies  of  the  levels  associated  with  the  ion  limit  ^PqI,  the  number 
is  further  reduced  to  five. 

Shortley  and  Fried  derived  approximate  theoretical  expressions 
for  the  energies  of  the  levels  of  the  p^d  and  p-*f  configurations  for  pure 
pair  coupling  in  terms  of  the  energy  parameters,  or  Slater  integrals, 
under  the  limiting  assumptions  that  only  the  leading  term  in  the  electro¬ 
static  interaction  and  the  splitting  factor  of  the  parent  configuration 
need  be  considered. 

Accordingly,  instead  of  the  usual  six  significant  energy  parameters, 
only  three  were  left,  namely: 


4 

NAVWEPS  Report  7190,  Quarterly  Report:  Foundational  Research 
Projects,  April-June  1961.  p.  11  (1961). 

5 

G.  H.  Shortley  and  B.  Fried,  Phys.  Rev.,  Vol.  54,  p.  749  (1938). 

6G.  Racah.  Phys.  Rev.,  Vol.  61,  p.  537  (1942). 

7 

NAVWEPS  Report  72  19,  Quarterly  Report:  Foundational  Research 
Projects,  January-March  1962,  p.  1  1  (1962). 
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F0.  giving  the  location  of  the  group  of  levels 

F 2<  the  electrostatic  interaction  integral,  leading  to  the  distribu¬ 
tion  of  the  double  levels 

4p,  the  splitting  factor**  of  the  parent  ion,  permitting  evaluation 
of  the  separation  of  the  two  level  groups 

The  tabulation  of  the  quantum  designations  of  the  levels  arising 
from  the  p^f  configuration  in  pair  coupling,  including  the  appropriate 
vector  summations,  together  with  the  energy  formulas  derived  by 
Shortley  and  Fried,  is  displayed  in  Table  1.  Undoubtedly  owing  to  the 
lack  of  illustrative  data,  Shortley  and  Fried  did  not  work  out  the  deri¬ 
vation  for  p^g.  In  fact,  at  the  time  they  had  little  to  work  on  to  verify 
their  treatment  of  the  p^f  configuration.  The  following  prophetic  state¬ 
ment,  interesting  in  retrospect  after  twenty-five  years,  is  quoted:  "The 
only  observations  of  p^f  are  in  argon  where  there  are  sketchy  data  for 
a  number  of  these  configurations.  At  most,  five  of  the  twelve  levels 
are  observed  and  not  even  the  J- values  of  these  are  certain.  It  may  be 
that  the  above  formulas  can  be  of  use  in  interpreting  these  observations. 
Perhaps  the  five  observed  levels,  four  lower  and  one  upper,  are  at  the 
positions  of  the  five  collapsed  levels  given  by  these  formulas." 

Atomic  Energy  Levels,  compiled  by  Moore,  ^  lists  both  experimen¬ 
tal  and  calculated  energy  levels  from  the  p^f  configurations.  The  latter' 
are  from  an  unpublished  manuscript  by  Edlen.  In  1956,  Eriksson*® 
published  a  paper  entitled  "Coupling  of  Electrons  with  High  Orbital 
Angular  Momentum,  Illustrated  by  2p  nf  and  2p  ng  in  N  II."  A  set  of 
formulas  for  the  pg  configuration  was  derived  by  a  treatment  essen¬ 
tially  identical  with  that  employed  by  Shortley  and  Fried.  Eriksson  had 
included  the  effect  of  interaction  between  levels  converging  to  different 
ion  limits.  This  included  an  extra  term  in  each  of  two  of  the  energy 
equations.  By  leaving  off  the  extra  term — something  later  justified  by 
treatment  of  the  data  —  and  changing  the  algebraic  signs  of  the  coeffi¬ 
cients  of  F2  and  £p  in  accordance  with  the  inversion  of  levels  associ¬ 
ated  with  complementary  equivalent- electron  configurations,  the  Eriks¬ 
son  expressions  convert  to  a  form  equivalent  to  those  derived  by  Short¬ 
ley  and  Fried.  Table  2  is  developed  for  the  p^g  configuration  in  the 
same  form  as  Table  l  for  the  p^f  configuration. 

A  detailed  treatment  of  the  derivation  of  the  energy  formulas  is 
not  included  here.  The  referenced  articles  may  be  consulted.  For  a 
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S.  Goudsmit  and  C.  J.  Humphreys,  Phys.  Rev.,  Vol.  31,  p.  960  (1928). 
9 

C.  E.  Moore,  Atomic  Energy  Levels,  Natl.  Bur.  Standards  Cir¬ 
cular  467  (1949). 

10K.  B.  S.  Eriksson,  Phys.  Rev.,  Vol.  102,  p.  102  (1956). 
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TABLE  1.  Shortley- Fried  Approximation  for  p  f 


(Taken 

from  Phya.  Rev., 

Vol.  54, 

p.  749,  1938) 
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5.4 
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7/2 
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Fo+S> 

1/2 

3 
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TABLE  2. 

5 

Approximation  for  p  g,  Based  on  Shortley- 
Fried  Treatment 

•*1 

1 

K 

J 

Energy 

3/2 

4 

11/2 

6,5 

F»-  l/2tp-28F2 

3/2 

4 

9/2 

5,4 
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4 
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4 
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Fo  +  tP 

3/2 

4 

7/2 

4,3 
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3/2 

4 

5/2 

3,2 
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more  extensive  discussion,  the  standard  work  of  Condon  and  Shortley^ 
is  suggested. 

Examination  of  the  energy  equations  shows  that  only  two  unknowns 
need  to  be  determined  for  a  solution  since  £p  is  related  to  the  separa¬ 
tion  of  the  ion  doublet  by  the  relation  Ap^  =  3/2£p.  Correct  identifi¬ 
cation  of  any  two  transitions  involving  levels  whose  energy  description 
includes  F^  permits  evaluation  of  the  parameters  and  derivation  of  the 
remaining  levels,  which  can  in  turn  be  checked  against  appropriate 
transitions.  Certain  clues  to  the  identification  of  the  selected  transi¬ 
tions  are  fairly  obvious.  Among  the  4f-5g  transitions,  4f[4j]^  - 
5g[  5j]  £  is  certain  to  be  the  most  intense.  Of  the  array  of  pairs  of  f- 
levels,  4f  [4|  4  and  4f[3|]4  3  have  the  smallest  pair  separations; 

in  fact,  neither  has  ever  been  resolved  experimentally.  This  condition 
leads  to  the  expectation  that  the  features  resulting  from  the  double  tran¬ 
sitions 


« l4i |  5  4  -  5g[5il-  5  and 

4lUtl4,3-5g[4||-4 

will  be  appreciably  sharper  than  the  others.  Finally,  based  on  the 
sequence  of  levels  in  each  category,  a  feature  common  to  the  spectra 
of  all  the  noble  gases,  the  most  intense  transitions  that  form  the  diago¬ 
nal  of  the  supermultiplet  come  in  the  following  order  of  increasing  wave 
number : 

«l3tl4jJ-5g[4ii;i4 

4tUil3,2-5g[3il;3 

4f[4tl5>4-5g[5il;i5 

4fHJi2il-5gU4i;t2 

These  characteristics  are  all  illustrated  by  Figure  2,  a  record  of  Ar  I 
in  the  4-micron  region  under  the  best  available  conditions  of  resolution. 
The  entire  procedure  becomes  an  exceedingly  simple  and  powerful  clas¬ 
sification  tool  in  instances  where  pair  coupling  applies. 


E.  U.  Condon  and  G.  H.  Shortley,  Theory  of  Atomic  Spectra, 
Cambridge  University  Press  (1951). 
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FIGURE  2.  Section  of  a  Recording  of  the  First  Spectrum  of  Argon  in  the 
4- Micron  Region,  with  Identifications  of  the  4f-5g  Transitions 


For  various  degrees  of  intermediate  coupling  where  the  other  Slater 
parameters,  namely,  the  exchange  integrals  and  the  spin  coupling  of  the 
valence  electron,  are  significant,  similar  treatments  are  useful  but  are 
much  more  involved,  requiring  iterated  approximations  to  arrive  at  the 
best  fit  between  the  calculated  eigenvalues  of  the  Hamiltonian  and  the 
experimental  energy  levels  of  the  atom. 

Although  this  report  is  primarily  concerned  with  the  evaluation  of 
levels  originating  in  the  p^g  configuration,  a  few  comments  are  inter¬ 
jected  regarding  the  character  of  the  vector  coupling  associated  with 
noble-gas  configurations  involving  smaller  f  -values  of  the  outer  elec¬ 
tron,  specifically  the  configuration  p^f.  In  all  the  noble  gases  some 
degree  of  intermediate  coupling  is  in  evidence  for  all  configurations  of 
greater  excitation  than  the  normal  state,  which  is  associated  with  s^p^. 

The  development  of  pair  coupling  proceeds  much  more  rapidly  for 
neon  than  for  the  other  noble  gases  and  is  essentially  complete  with  the 
configuration  p^f.  The  first  detailed  and  complete  description  of  the 
4f-level8  was  based  on  the  observation  of  the  3d-4f  transitions  by 
Humphreys  and  Kostkowski.  ^  Eriksson  10  used  this  material  to  illus¬ 
trate  his  theoretical  treatment  but  presented  it  only  in  graphic  form. 

We  are  therefore  inserting  a  tabulation  of  the  experimental  4f-levels 
and  the  values  calculated  on  the  basis  of  the  pair -coupling  approxima¬ 
tion,  according  to  the  equations  in  Table  1.  This  material  is  shown  in 
Table  3.  Also  listed  are  the  A.  E.  L.  values,  including  unobserved 
values  calculated  by  Edlen.  The  only  significant  difference  between  the 
observed  values  and  those  calculated  on  the  basis  of  the  pair-coupling 
hypothesis  is  that  two  coalesced  pairs  associated  with  the  limit 

are  observed,  instead  of  a  complete  collapse  into  a  single  level  of  four 
components . 

It  is  noted  also,  with  reference  to  Table  3  and  the  level  tabulations 
for  Ar  I  and  Kr  I  that  follow,  that  the  calculated  values  depend  some¬ 
what  on  which  experimental  pair  of  levels  is  introduced  into  the  calcu¬ 
lations.  A  final  publication  is  expected  to  include  a  description  based 
on  more  measurements  and  adoption  of  level  values  derived  from  an 
average  of  results  obtained  from  alternative  sets  of  starting  data.  A 
clue  to  the  degree  of  pair  coupling  is  to  be  found  from  examination  of 
the  four  collapsed  pairs  converging  to  the  limit  ^P}I.  Because  of  the 
values  of  the  coefficients  of  F2  for  the  various  pairs  of  levels  in  pure 
pair  coupling,  when  the  four  energy  values  are  arranged  in  a  monotonic 
sequence  the  upper  and  lower  pairs  have  exactly  the  same  interval 
separating  them. 


C.  J.  Humphreys  and  H.  J.  Kostkowski,  J.  Research  Natl.  Bur- 
Standards,  Vol.  49,  p.  73  (1952). 
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‘Based  on  Humphreys  -  Kostkow  ski  values. 


This  condition  is  realized  to  well  within  the  experimental  error  for 
these  four  levels  of  Ne  1.  Increasing  departure  is  noted  for  the  other 
noble  gases.  In  Xe  I  the  4f-levels  fit  the  pair-coupling  relationship 
rather  poorly.  Another  indication  is  the  separation  of  the  pair  having 
the  double  J-value  1,  2.  Unresolvable  in  Ne  I,  it  is  interferometrically 
separable  in  Ar  I  with  an  interval  of  0.40  cm-*,  it  is  easily  resolvable 
in  Kr  I  with  an  interval  of  0.88  cm'',  but  in  Xe  I  it  shows  a  relatively 
wide  separation  of  9.67  cm''.  Similarly,  the  4f-level  pair  having  the 
double  J-value  2,  3,  unresolvable  in  Ne  I  and  Ar  I,  is  resolved  in  favor¬ 
able  instances  in  Kr  I  and  separated  by  3.0  cm"  *  in  Xe  I.  The  4f-level 
with  J-values  4,  5  is  resolvable  only  in  Xe  I. 

In  order  to  compare  the  characteristic  features  of  the  multiple 
group  of  lines  associated  with  the  4f-5g  transitions,  a  set  of  sections 
of  records  for  the  various  noble-gas  spectra,  aligned  as  to  wavelength, 
has  been  assembled  as  Figure  3.  Particularly  to  be  noted  are  the  shift 
of  the  multiple  group  to  shorter  wavelengths  with  increasing  atomic 
number,  and  the  transition  from  a  broadening  of  unresolved  features  to 
partial  and  finally  complete  resolution,  also  with  increasing  wave  num¬ 
ber.  With  reference  to  transitions  between  levels  converging  to  the 
ion  limit,  a  partially  resolved  feature  is  noticeable  only  in  Ar  I. 

Of  the  two  pairs  of  4f  -levels  in  this  spectrum,  one  is  much  closer  than 
the  other,  0.35  cm''  for  the  f*  3  3  pair  as  against  0.05  cm"  *  for  the 
f  3,4  pair.  This  accounts  for  the  double  feature  with  one  sharp  and  one 
diffuse  component.  In  the  instance  of  Kr  I,  the  5g*  complex  is  above 
the  first  ion  limit,  and  in  Xe  I,  both  the  4f’  -  and  5g'  -levels  are  above 
this  limit.  In  both  instances,  although  the  structure  should  be  resolv¬ 
able,  only  one  feature  is  observed  that  is  probably  associated  with  the 
most  intense  transition. 

Taking  up  now  the  consideration  of  the  5g-levels  determined  from 
the  4f-5g  transitions,  observation  of  which  has  supplied  the  experimen¬ 
tal  material  for  this  study,  it  is  found  that  these  levels  in  Ne  I,  Ar  I, 
and  Kr  I  all  conform  to  the  pair -coupling  approximation  within  the  limits 
of  observational  error. 

The  corresponding  Xe  I  transitions  have  been  observed  but  it  has 
not  been  found  possible  to  sort  out  the  levels  by  a  similar  treatment. 
Evidently  the  departure  from  pair  coupling  is  too  great  in  this  instance 
so  that  a  more  detailed  treatment  involving  additional  parameters  will 
be  required.  This  will  be  reserved  for  a  future  report. 

The  results  of  the  determinations  of  the  5g-levels  of  Ne  I,  Ar  I, 
and  Kr  I  are  now  presented  in  order  by  means  of  arrays  shown  as 
Tables  4,  5,  and  6.  In  each  instance  the  4f-levels  are  shown  in  the 
left-hand  column  in  the  same  order.  The  g-levels  are  arranged  mono- 
tonically  at  the  top  of  the  page.  The  arrangement  of  the  set  of  double 
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FIGURE  3.  Superposed  Sections  of  Recordings  of  the  Spectra  of  Neon, 
Argon,  Krypton,  and  Xenon,  with  Correlation  of  Analogous  Features 
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f-level  of  largest  J  used  to  calculate  v  in  transitions. 


J-values  with  respect  to  the  monotonic  sequences  is  uniform  for  all  the 
noble  gases  and  the  pattern  is  maintained  for  both  the  4f  and  5g  arrays. 

A  summary  compilation  of  the  values  of  the  energy  parameters 
for  the  5g  configurations  for  the  elements  neon,  argon,  and  krypton, 
as  derived  on  the  basis  of  the  selected  starting  data,  appears  in  Table  7. 
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TABLE  7.  Energy  Parameters  for  the  p  g 
Configuration 


Element 

Cp.  cm"1 

F q,  cm'  1 

F 2»  cm"  1 

Neon 

520.33 

169801.93 

0.05192 

Argon 

954.70 

123193.87 

0.13363 

Krypton 

3580.67 

110310.16 

0.17701 

In  some  noted  instances,  other  Ritzian  combinations  appear  on  the 
records.  Some  of  these  may  be  seen  in  the  various  sections  of  Figure 
3.  The  availability  of  such  transitions  is  advantageous  for  data  reduc¬ 
tion  because  the  wavelengths  associated  with  them  are  precisely  cal¬ 
culable.  In  moat  instances  the  wavelengths  were  measured  by  intro¬ 
ducing  a  superposed  second-order  argon  spectrum. 

Finally,  some  observations  have  been  made  on  He  I  by  using  a 
source  tube  containing  a  helium-neon  mixture.  Figure  4  displays  a 
record  from  this  source  superposed  on  one  from  a  pure  neon  source. 
Three  features  attributed  to  He  I  are  identified  as: 


X,  air,  A 

cr,  vac.,  cm"' 

Transition 

40364.2 

2476.77 

4d3D 

-  5f3F 

40395.6 

2474.84 

4dJD 

-  5f*F 

4f3F* 

-  5g3< 

40477.9 

2469.81 

4fJF* 

-  5g  i< 

The  wavelengths  of  the  first  two  listed  helium  lines,  particularly 
the  singlet  transition,  differ  somewhat  from  those  calculated  from  level 
values  as  published  in  Atomic  Energy  Levels.  It  is  expected  that  some 
revision  of  level  values,  as  determined  by  further  measurements,  will 
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FIGURE  4.  Section  of  a  Recording  of  the  Spectrum  of  a  Helium-Neon 
Mixture  Compared  with  a  Recording  of  the  Spectrum  of  Neon  Covering 
the  Same  Wavelength  Interval  Near  4  Microns 


be  required.  The  interpretation  of  the  third  feature  is  not  verifiable  by 
intercombinations  and  possible  ambiguity  is  indicated.  The  validity  of 
the  interpretation  of  the  triplet  combination  seems  unquestionable  owing 
to  the  nearly  hydrogenic  character  of  the  high  levels.  It  is  noted  that 
the  Brackett  line  of  hydrogen  has  a  wavelength  of  40511  A. 

The  authors  are  indebted  to  K.  L.  Andrew  of  Purdue  University  for 
suggesting  the  use  of  the  Shortley-Fried  approximation  in  the  interpre¬ 
tation  of  the  experimental  data,  and  to  R.  D.  Cowan  of  the  University  of 
California,  Los  Alamos  Scientific  Laboratory,  for  valuable  discussions. 
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NONAQUEOUS  ELECTROCHEMISTRY 


THE  ELECTROREDUCTION  OF  NITROBENZENE 
IN  LIQUID  AMMONIA  SOLUTIONS 

by 

W.  S.  Harris 

For  some  time  it  has  been  known  qualitatively  that  the  m-dinitro- 
benzene  (abbreviated  here  as  m-DNB)  cathode,  when  used  in  current- 
producing  cells  and  batteries,  performs  better  in  an  acid  (NH4SCN) 
than  in  a  neutral  (KSCN)  environment.  This  phenomenon  was  investi¬ 
gated  as  a  part  of  NOLC's  study  of  the  reduction  mechanism  of  aromatic 
nitro  compounds  in  liquid  ammonia  solutions.  The  results  from  this 
investigation  were  utilized  to  recalculate  the  maximum  energy  to  be 
expected  from  the  electrochemical  cell  Li/electrolyte  in  NH^/m-DNB. 


EXPERIMENTAL 

The  apparatus  and  method  are  the  same  as  reported  in  previous 
quarterly  reports  ^  and  in  NAVWEPS  7241. ^ 


RESULTS 

rn-Dinitrobenzene  behaves  quite  differently  when  reduced  in  acid 
liquid  ammonia  solutions  than  when  reduced  in  neutral  liquid  ammonia 
solutions.  This  behavior  is  shown  in  Figure  la.  At  low  rates,  in  acid 
solutions  the  useful  yield  is  the  expected  8  electrons  per  molecule;  in 
neutral  solutions  the  useful  yield  is  about  2  electrons.  The  useful  yield 


NAVWEPS  7214,  Quarterly  Report:  Foundational  Research  Pro¬ 
jects,  October-December  1961,  pp.  29-35  (March  1962);  and  NAVWEPS 
7229.  Quarterly  Report:  Foundational  Research  Projects,  April-June 
1962,  pp.  33-39  (August  1962). 

2 

NAVWEPS  7241,  The  Electroreduction  of  Aromatic  Nitro  Com- 
pounds  in  Nonaqueous  Solutions.  Part  I;  The  Reduction  of  Nitrobenzenes 
in  Acid  Liquid  Ammonia  Solutions,  by  W.  S.  Harris  and  G.  B.  Matson 
(December  1962). 
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FIGURE  la.  Low-Rate  Reduction  of  m-DNB  in  Liquid  Ammonia  Solutions 


ATURE 


FIGURE  lb.  Low-Rate  Reduction  of  m-DNB  in  Liquid  Ammoni?  Solutions 


is  the  number  of  electrons  used  until  a  0. 4-volt-*  drop  in  cathode  poten¬ 
tial  occurs. 

The  addition  of  strong  acids  (in  ammonia),  such  as  NH4SCN,  to  the 
neutral  solution  increases  the  utilization  of  electrons  approximately  in 
proportion  to  the  amount  of  acid  added.  When  a  weaker  acid  is  added, 
a  much  greater  amount  must  be  used  to  give  the  same  effect.  This  is 
illustrated  in  Figure  lb.  Cells  651  and  652  gave  similar  results  but 
contained  different  relative  amounts  of  added  acid.  Cell  652  used  the 
strong  acid  NH4SCN  with  an  acid:m-DNB  mole  ratio  of  4.  Cell  651 
used  the  weak  acid  H2O  at  a  ratio  of  100. 

The  addition  of  sufficient  strong  acid  to  the  neutral  solution  to  bring 
the  acid:m-DNB  mole  ratio  to  10  results  in  a  discharge  curve  that  is 
virtually  the  same  as  in  acid  solution,  with  the  exception  of  a  slight 
drop  in  voltage. 


CONCLUSIONS 

From  these  investigations,  several  conclusions  are  drawn  for  well- 
stirred  solutions  at  very  low  current  densities: 

(1)  m-DNB  is  incompletely  reduced  in  neutral  solutions. 

(2)  Acid  is  needed  for  the  complete  8-electron  reduction  of  m-DNB. 

(3)  The  concentration  of  acid  appears  to  be  of  minor  importance. 

It  is  postulated  that  these  conclusions  will  have  the  following  effects 
upon  hardware  type  cells: 

(1)  rn-DNB  will  be  utilized  best  in  cells  having  a  NH4+:m-DNB 
mole  ratio  of  8  or  higher. 

(2)  Cells  containing  smaller  amounts  of  NH4+  than  this  ratio  may 
be  expected  to  have  variable  lives  and  voltages  that  depend 
upon  the  amount  of  NH4+  present  in  the  cathode. 

(3)  Because  NH4+  is  an  essential  part  of  the  reaction,  it  must  be 
included  in  any  weight  calculation  that  includes  the  m-DNB 
cathode. 


This  value  is  arbitrary  and  cannot  be  related  to  battery  voltage 
regulation. 
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CALCULATIONS 


Because  H+  participates  in  and  is  necessary  for  the  reaction 
m-<f>(N02)2 - ►m-((>(NHOH)2 

it  must  be  included  in  any  energy  density  calculation  for  cells  utilizing 
rn-DNB  as  a  cathode  material.  Energy  density  calculations  for  the  sys 
tern  Li/electrolyte  in  NH3/m-DNB,  based  upon  various  assumptions, 
are  given  below. 

Calculation  No.  1 
Assumptions : 

(1)  Need  for  is  disregarded. 

(2)  100%  efficiencies. 

(3)  Working  voltage  is  same  as  E°  (2.6  volts,  estimated). 

Cell  Reaction: 

8  Li  +  rn-DNB  — — — ►  8  Li+  +  other  products 
Watt-hours/pound^  =1100 

Calculation  No.  2 
As  sumptions : 

(1)  Source  of  H+  is  NH};  i.  e.,  NH3  - r  H+  +  NHf 

(H^O  would  give  a  similar  result.  ) 

(2)  100%  efficiencies. 

(3)  W  orking  voltage  is  same  as  E°. 

Cell  Reaction: 

8  Li  +  m-DNB  +  8  NH3  — — — ►  8  UNH2  +  m-<)>(NHOH)2  +  2  H20 
Watt-hours/pound  =  705 

Calculation  No.  3 
As  sumptions : 

(1)  Source  of  H+  is  NH4+  from  NH4F.5 

(2)  100%  efficiencies. 

(3)  Working  voltage  is  same  as  E°. 


One  watt-hour /pound  =  7. 04  joules/gram. 
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NH4F  was  arbitrarily  taken  as  one  of  lightest  H+  sources.  It  is, 
however,  insoluble  in  NH3. 


Cell  Reaction: 

8  Li  +  m-DNB  +  8  NH4F  - 6  ►  8  LiF  +  m-<j>(NHOH)2 

+  8NH3  +  2H20 

Watt-hours/pound  =  482 

Calculation  No.  4 
Assumptions: 

(1)  Source  of  H+  is  the  NH4+  from  NH4NO3. 

(2)  100%  efficiencies. 

(3)  Working  voltage  is  same  as  E°. 

Cell  Reaction: 

8  U  +  m-DNB  +  8  NH4NO3  +  m-<j>(NHOH)2  +  8  NH3  +  2  H20 
Watt-hour  s/pound  =  295 

Calculation  No.  5 
Assumptions : 

(1)  Source  of  NH4+  is  NH4NO3. 

(2)  50%  efficiencies. 

(3)  Corrosion  reaction  at  anode  causes  its  inefficiency. 

Li  +  NH4+ - ►  Li+  +  1/2  H2  +  NH3 

(4)  Lack  of  reaction  causes  inefficiency  at  cathode. 

(5)  Working  voltage  is  2.2  volts  at  100  ma/cm^. 

Cell  Reactions : 

8  Li  +  m-DNB  +  8  NH4NO3  - ►  8  UNO3  +  m-4»(NHOH)2 

+  8NH3  +  2H20 

8  Li  +  8  NH4NO3 - ►  8  LiN03  +  4  H2  +  8  NH3 

rn-DNB - ►rn-DNB  (i.e.,  no  reaction) 

Watt-hours/pound  =  122 

Calculation  No.  5  is  the  best  approximation  for  the  system  as  it  is 
now  known  in  research-hardware  type  cells  at  room  temperature. 
Assumption  (3)  is  based  upon  a  Li  anode  in  a  NH4NO3-NH3  environment. 
The  corrosion  would  be  much  less  in  a  neutral  solution. 


From  the  above  calculations,  it  is  seen  that  the  dinitrobenzene 
cathode  is  penalized  severely  because  of  its  need  for  acid.  To  achieve 
better  energy  it  is  necessary  to  find  (1)  strong  acids  with  lower  equiva- 


lent  weights  or  (2)  cathode  materials  that  do  not  require  acid  or  that 
can  use  weaker  acids  such  as  H2O. 


IONIC  MELT  ELECTROCHEMISTRY 
by 

R.  E.  Panzer 


Investigations  in  ionic  melts  reported  in  the  last  quarterly  report 
have  been  continued.  Problems  of  instrumentation  and  preparing  elect¬ 
rolyte  pills  were  solved  and  all  necessary  cell  components,  i.e.,  anodes, 
anolyte  pills,  catholyte  pills,  cathodes,  and  insulators,  were  prepared. 

Fifty-four  factorial  cell  experiments  were  run  at  the  rate  of  18 
cells  per  day.  Each  daily  set  was  operated  at  one  of  three  tempera¬ 
tures:  200,  250,  or  300*C.  These  temperatures  differ  from  the  origi¬ 
nal  plan  of  175,  250,  or  325*C,  the  change  being  necessitated  by  prob¬ 
lems  encountered  in  preliminary  cell  tests.  Statistical  evaluation  and 
analysis  of  the  data  from  these  cell  tests  are  in  process. 

A  literature  search  for  prospective  low  melting  electrolytes  is 
continuing.  Several  mixtures  of  thiocyanates  and  other  salts  appeared 
quite  promising  in  preliminary  evaluation  of  melting  point,  conduc¬ 
tivity,  stability,  and  unreactivity  with  cell  components.  Results  of 
these  studies  will  be  reported  later. 

In  conjunction  with  studies  of  electrochemical  polarization  in  prac¬ 
tical  cell  systems,  efforts  are  being  made  to  apply  modified  techniques 
of  voltammetry  at  constant  current  (chronopotentiometry) .  Problems 
under  study  are  reproducible  electrode  areas  and  their  positioning  in 
the  cell,  and  concentration  limits  in  ionic  melt  voltaic  cells.  Circuitry 
has  been  developed  and  tested  by  using  short  time  (<  1  second)  electrol¬ 
ysis  with  oscilloscope  recording.  Additionally,  longer  electrolysis 
times  of  up  to  5  seconds  have  been  used  with  a  fast  strip  chart  recorder. 
Presently,  it  is  not  known  whether  this  technique  will  yield  interpre¬ 
table  results  in  the  practical  ionic  melt  cells. 
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NONLINEAR  TRANSMISSION  LINES 


THE  SEMIGENERALIZATION  OF  THE 
TRANSMISSION  LINE  EOUATIONS 

by 

J.  R.  Alday 


INTRODUCTION 

In  previous  studies,  attempts  have  been  made  to  use  the  general 
linear  form  of  the  transmission  line  equations  to  describe  the  actions 
of  a  nonlinear  transmission  line.  A  lumped  constant  approach  has  also 
been  attempted.  Both  of  these  approaches,  although  quite  fruitful  in 
some  instances,  are  insufficient  for  fully  describing  the  nonlinear  line. 

The  analysis  presented  here  is  an  attempt  to  generalize  further 
the  usual  linear  transmission  line  equations  to  include  the  case  of 
active  or  magnitude-sensitive  admittance  elements.  The  term  "semi¬ 
generalization"  has  been  coined  for  use  where  only  the  admittance  ele¬ 
ments  are  considered  active. 


ANALYSIS 

An  equivalent  circuit  of  infinitesimal  length  (Figure  1)  is  used  here 
to  develop  the  nonlinear  transmission  line  equations  in  the  same  man¬ 
ner  as  it  was  used  for  the  linear  transmission  line  equations. 


Z 

— WV 


y(v(I)]=  y[i  (V)] 


dz- 


FIGURE  1.  Equivalent  Circuit 
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Let  Z  be  the  complex  series  impedance  of  the  infinitesimal  section 
of  the  transmission  line  and  let  Y[V(I)]  or  Y[I(V)]  be  a  functional  repre¬ 
sentation  for  the  complex  shunt  admittance  in  I  or  V  respectively.  The 
terms  V(I)  and  I(V)  are  indicative,  respectively,  of  the  voltage  across 
and  the  current  through  the  admittance  element. 

From  Ohm's  Law,  the  following  equation  may  be  written: 

V'IZT  =  I|Z+YTf(I>lj  "> 

where  V  is  the  instantaneous  input  voltage  signal,  Zf  is  the  total  imped¬ 
ance  of  the  equivalent  circuit,  and  I  is  the  instantaneous  current  drawn 
by  the  circuit. 

If  V  is  allowed  to  change  to  V  +  AV,  there  will  be  a  corresponding 
change  in  I;  that  is 

V  ♦  AV  .  (I  ♦  AI>Jz  ♦  jj—V—jj  ,2) 

Expanding  equation  (2)  and  subtracting  equation  (1)  from  it  gives 
AV  1  1  .  AT7  i 

Y  [  V(I  +  AI)]  '  Y[V(I)]  a  Y  [  A  (I  +  AI)] 

Collecting  terms  and  dividing  both  sides  by  AI  yields 

AV  _  <Y[V(I)j  -  Y [ V(I  +  AI)]  j  1 

AI  |AI-Y[V(I  +  AI)]  Y[V(I)]  (  Y[V(I+AI)] 

Observing  that 
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Y  [  V(I) J  -  Y  [  V(I  +  AI)]  =  -  AY  [  V(I)] 
and  passing  to  the  limit  results  in 

lim  *1  s  il  -i  iXLXffll _ ! _ *z+  _ 

Al-0  AI  dI  dI  Y2[V(I|]  Ytv«I>] 

Integrating  *  with  respect  to  I  yields 

I 


V  = 


Y[V(I)] 


+  ZI  + 


/ 


dl 


Y[v(i)] 


^The  integration  at  the  first  term  on  the  right-hand  side  is  actually 
with  respect  to  Y[V(I)j  . 


S' 


Now,  inasmuch  as  the  change  in  voltage  per  unit  length  of  the  infini¬ 
tesimal  section  of  line  is  equal  to  the  voltage  drops  around  the  circuit, 

8V (I)  I  _  8V  ?r  f  dl 
8z  ~  Y[V(I)]  8z  J  Y[V(I)] 

Also,  since  the  applied  signal  V  is  not  a  function  of  z,  its  derivative 
vanishes  and 

TT1=-/lztYTW)l|dI  ,3) 

Equation  (3)  is  the  semigeneralized  form  of  the  voltage  transmission 
line  equation. 


In  a  very  similar  manner  (the  details  of  which  are  not  included  here), 
the  semigeneralized  current  transmission  line  equation  is  found  to  be 


dV 


1 

Y[  I(V)j 


The  usual  equations 


and 


8V 

8z 


=  -ZI 


8I_ 

8z 


=  -YV 


(4) 


(3a) 

(4a) 


where  Z  and  Y  are  complex,  can  be  similarly  derived  if  it  is  assumed 
that  Y  is  invariant  with  voltage. 

For  further  generalization,  assume  that  the  time -dependent  volt¬ 
ages  and  currents  are  Fourier  in  nature,2  that  is  to  say,  they  take  on 
the  respective  forms 


00 

1=  E  Inexp(i[nu,t  +  <>n(t)]}  0<t<-^  (5) 


and 


This  assumption  is  not  only  academically  important.  Its  physical 
interpretation  is  at  least  this  complex. 
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(6) 


00 

V  =  2  V  exp{i[nwt  +  0  (t)l } 
,  n  n 

n=l 


0  <  t  < 


Ztt 

w 


where 

n  =  If  Z ,  3,  •  •  • 

In  =  maximum  current  magnitudes  of  the  fundamental  signal  and  its 
harmonics 

Vn  =  maximum  voltage  magnitudes  of  the  fundamental  signal  and  its 
harmonics 

4>n(t)  =  complex  current  phase  functions 
0n(t)  =  complex  voltage  phase  functions 

Direct  substitution  of  equations  (5)  and  (6)  into  equations  (3)  and  (4), 
with  the  change  in  variables  Y[V{I)]  =  Y[I(V)]  =  Y(t),  results  in  the 
semigeneralized  equations 


av(t) 

9z 


f  00  r  84>(tn 

=  -  /  Z(t)  £  Ini  nw  +  — — —  I  exp  {i[nwt  +  <t>n(t)] }  dt 
n=l  L  ■* 


(7) 


and 


9I(t)  _ 
9z 


/.  i  00  T  90  (t)  "I 

=-/^,vK] 


exp{i[nu)t  +  0  (t) ] }  dt 


(8) 


where 


Z(t)  =  Z  +  Y[V(I)]  =  z  +  y(t) 


CONCLUSIONS 

The  transmission  line  equations  for  a  transmission  line  with  active 
shunt  admittances  have  been  derived.  It  is  felt  that  these  equations, 
although  more  complex  than  those  that  do  not  consider  the  admittance 
elements,  are  representative  of  the  properties  of  a  nonlinear  trans¬ 
mission  line.  Complete  utilization  of  these  equations  is  dependent  upon 
their  simultaneous  solution. 
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SEMICONDUCTOR  MICROWAVE  MICROELECTRONICS 


SUBSTRATES.  CIRCUITRY.  AND  SEMICONDUCTORS 

by 

R.  W.  Yancey 


The  last  quarterly  report  contained  a  brief  discussion  of  some 
aspects  of  microwave  semiconductor  microelectronics.  In  continuation 
of  this  study  to  develop  and  determine  properties  of  microwave  func¬ 
tional  circuits  employing  semiconductors,  tunnel  diodes  were  chosen 
as  the  major  semiconductor  device  to  be  considered.  The  properties 
and  advantages  of  tunnel  diodes  have  been  covered  in  most  electronic 
publications  and  need  not  be  repeated  here.  These  devices  are  capable 
of  operation  beyond  100  Gc.l  but  commercial  units  are  limited  to  30  Gc. 
This  30-Gc  limitation  precludes  the  use  of  tunnel  diodes  in  millimeter 
wave  equipment  which  has  a  number  of  ordnance  as  well  as  other  appli¬ 
cations. 

Consideration  was  given  to  taking  advantage  of  the  tunnel  diodes' 
capabilities  by  fabricating  the  tunnel  diodes  directly  into  microwave 
structures  to  increase  the  frequency  of  operation  into  the  millimeter 
wave  region.  A  brief  search  revealed  that  tunnel  diodes  are  made  by 
the  alloying  process,  which  is  simple  and  straightforward  if  an  alloy¬ 
ing  furnace  is  available. 

A  trip  to  the  Harry  Diamond  Laboratories  (which  was  then  the 
Diamond  Ordnance  Fuze  Laboratory)  resulted  in  the  acquisition  of  suf¬ 
ficient  information  to  build  an  alloying  furnace  and  to  perform  the  alloy¬ 
ing  process.  (It  should  be  mentioned  that  these  first  costs  were  borne 
by  Code  71  during  one  of  its  programs  that  was  later  terminated.  The 
work  was  then  continued  under  Foundational  Research  funds. ) 

The  decision  was  made  to  go  ahead  in  this  direction  and  to  fabri¬ 
cate  the  diodes  into  microwave  equipment.  From  the  information  and 
sample  material  obtained  from  DOFL,  the  furnace  was  constructed  as 
follows . 


C.  A.  Burrus,  "Gallium  Arsenide  Esaki  Diodes  for  High-Fre¬ 
quency  Applications,"  J,  Appl.  Phys.,  Vol.  32,  pp.  1031-1036  (1961). 
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The  furnace  (Figure  1)  was  made  of  stainless  steel  and  glass,  with 
a  base  7  inches  in  diameter  and  approximately  4  inches  high.  The  heat¬ 
ing  element  is  a  filament  made  from  a  0.010-inch  nickel  sheet  5/8  of 
an  inch  wide.  The  edges  are  bent  up  l/l6  inch  at  a  90-degree  angle  to 
give  rigidity.  The  two  leads  on  the  left  of  the  figure  are  external  con¬ 
nections  for  thermocouple  leads.  The  furnace  is  evacuated  through 
the  center  hole  in  the  base  and  the  reduction  atmosphere  is  admitted 
through  the  same  hole.  This  is  accomplished  by  an  external  glass 
"T."  The  filament  current  is  supplied  via  its  supporting  posts  which 
go  through  Teflon  insulators  to  the  base. 

The  closed  furnace,  with  an  over-all  height  of  6  inches,  is  shown 
in  Figure  2.  The  handles  exert  pressure  on  the  glass  flange  to  effect 
a  seal  from  an  O-ring  between  the  base  and  bell  jar.  The  furnace 
proper  is  now  complete. 

The  filament  temperature-control  unit  is  complete  with  the  exception 
of  the  timer  needed  to  give  completely  automatic  alloying  cycles.  This 
control  unit  is  capable  of  supplying  4  volts  at  75  amp  continuously,  and 
125  amp  intermittently.  The  anticipated  heater  current  is  30  to  50 
amp.  The  timer  will  be  adjustable  to  allow  30  to  50  seconds  for  pre¬ 
heating  and  2  to  10  seconds  for  alloying.  The  control  panel  is  shown 
in  Figure  3,  while  Figure  4  shows  a  rear  view  of  the  control  unit. 

Briefly,  the  alloying  process  consists  of  placing  the  doped  semi¬ 
conductor  wafer  on  a  heat-generating  element.  On  this  wafer  is  placed 
the  proper  alloying  material.  Both  are  then  placed  in  a  reducing  atmos¬ 
phere  (for  germanium,  gallium  arsenide,  or  gallium  antimonide,  hydro¬ 
gen  is  used)  and  preheated  to  550  or  600*C  for  approximately  45  sec¬ 
onds.  The  temperature  is  then  increased  to  approximately  700*C  for 
2  seconds,  and  the  unit  is  allowed  to  return  to  room  temperature.  The 
high  temperature  produces  the  alloying  which,  when  used  with  proper 
semiconductor  wafers  and  alloying  material,  yields  tunneling  junctions. 

Normally,  a  very  highly  doped  semiconductor  wafer  must  be  used 
to  produce  germanium  tunneling  junctions.  A  doping  level  on  the  order 
of  2  to  6  x  1019  atom^cm^  gives  a  resistivity  of  about  0.006  ohm-cm. 

The  wafer  thickness  may  vary  from  0.002  to  0.010  inch.  With  n-type 
germanium,  a  common  alloying  material  is  indium.  A  small  (0.0005- 
inch)  sphere  or  wafer  may  be  used. 

In  initial  tests,  4  germanium  tunnel  diodes  were  made  at  a  reduced 
pressure  (approximately  0.001  mm  Hg)  without  using  the  reduction  at¬ 
mosphere.  Since  no  provisions  were  made  for  mounting  these  units, 
they  are  at  best  just  samples. 
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FIGURE  1.  Alloying  Furnace 
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FIGURE  2.  Alloying  Furnace  with  Cover 
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FIGURE  4.  Back  View  of  Temperature -Control  Unit 


During  the  next  quarter,  the  furnace  control  will  be  completed, 
tested,  and  calibrated.  Full  operation  of  the  completed  system  is 
anticipated  by  the  next  reporting  period. 


SEMICONDUCTOR  PHYSICS 


ACCURATE  MEASUREMENTS  OF  A  MAGNETIC  FIELD 
BY  MEANS  OF  A  III-V  HALL  EFFECT  DEVICE 

by 

K.  I.  Crowder  and  H.  Piller 


INTRODUCTION 

A  current  traveling  through  a  material  in  the  presence  of  a  magnetic 
field  experiences  a  Lorentz  force  normal  to  the  magnetic  field  and  cur¬ 
rent  in  the  direction  J  x  B.  Because  the  material  constrains  the  elec¬ 
trons,  an  equal  and  opposite  electromotive  force  is  set  up,  2  =  Rif  x  T. 
This  electromotive  force  is  the  Hall  voltage.  Here  R  is  the  Hall  coef¬ 
ficient.  In  a  semiconductor  with  a  high  electron  concentration,  R  is 
inversely  proportional  to  this  concentration,  n,  and  therefore  R  =  l/ne, 
where  e  denotes  the  electronic  charge.  If  the  charge  carriers  are  posi¬ 
tive,  R  is  taken  to  be  positive  and  the  analysis  still  holds.  The  sense  of 
5  is  thus  determined  by  the  sign  of  the  carriers.  The  equation  for  R 
holds  in  III-V  compound  semiconductors  since  the  density  of  states  in 
the  conduction  band  is  low  and  the  contribution  of  electrons  to  the  con¬ 
duction  band  by  ionized  impurity  atoms  is  high  in  the  temperature  region 
where  the  Hall  generator  is  used. 

A  Hall  generator  has  many  advantages  over  other  methods  of  meas¬ 
uring  magnetic  field  strength.  In  the  presence  of  a  magnetic  field  for  a 
given  current,  it  generates  a  voltage  that  is  proportional  to  the  current, 
J,  and  also  the  magnetic  field,  B.  Its  physical  size  is  small;  therefore, 
it  samples  a  small  portion  cf  the  field.  The  output  voltage  is  determined 
by  the  average  field  strength  over  the  area  of  the  sample;  hence  the  field 
may  be  measured  in  small  gaps  and  small  areas. 


DISCUSSION 

The  voltage  output  is  influenced  by  temperature,  shunt  resistance, 
input  current,  and  magnetic  field.  The  temperature  dependence  of  the 
Hall  constant  R  is  equal  to  or  less  than  0.1%  per  degree  centigrade  in 
the  region  between  0  and  100*C.  By  using  manganin  wire  for  the  shunt 


67 


resistance  and  large  capacity  batteries  to  supply  the  input  current,  the 
accuracy  of  the  device  may  well  be  limited  only  by  the  temperature  sta¬ 
bility.  Error  due  to  this  instability  may  be  minimized  by  use  of  heat 
sinks  such  as  copper  tubes  with  water  from  a  temperature-controlled 
bath  circulating  through  them.  It  should  be  noted  that  orientation  in  the 
field  changes  the  apparent  thickness  of  the  device  by  a  factor  of  secant 
0,  where  0  is  measured  between  the  lines  of  magnetic  force  and  the  face 
of  the  Hall  generator.  The  voltage  is  proportional  to  cosine  0. 

A  Hall  generator  was  obtained  commercially  from  Instrument  Sys¬ 
tems  Corporation  (Siemens  Type  FA- 21)  and  calibrated  for  the  region 
between  2  and  10  kilogauss.  A  Leeds  and  Northrup  potentiometer, 

Model  K-3,  was  used  to  measure  the  output  voltage,  and  a  Harvey-Wells 
Corporation  nuclear  magnetic  resonance  precision  gaussmeter  system, 
Model  G502,  was  used  to  measure  the  field.  (A  schematic  diagram  of 
the  arrangement  is  shown  in  Figure  1.)  The  instrument  had  a  built-in 
crystal  calibrator  which  gave  resonance  lines  at  each  megacycle  fre¬ 
quency.  Because  the  field  could  be  determined  quite  accurately  (±0.01%) 
at  these  points,  they  were  chosen  for  the  calibration  of  the  Hall  genera¬ 
tor.  The  H*  and  the  Li^  probes  were  used. 

The  Hall  voltage  as  a  function  of  the  magnetic  field  deviates  from 
a  linear  behavior  in  a  uniform  and  predictable  manner  for  a  given  shunt 
resistance.  The  curve  has  an  S  shape  that  crosses  a  line  through  the 
origin  at  approximately  5  kilogauss  and  at  10  kilogauss  (see  Figure  2). 
This  curve  has  the  form 


V  =  aB  +  b  sin  (0.2  irB) 
h 

where  is  the  Hall  voltage  in  millivolts,  B  is  the  magnetic  field  in 
kilogauss,  and  a  and  b  are  constants.  The  sine  term  is  a  small  pertur¬ 
bation  on  the  linear  behavior  of  the  Hall  voltage  vs.  field;  b  will  not 
exceed  1%  of  the  value  of  the  Hall  voltage  at  10  kilogauss  if  the  proper 
shunt  resistance  is  chosen.  For  this  calibration 

a  =  4.661  millivolts/kilogaus  s 
b  =  3.867  millivolts 

The  best  accuracy  is  obtained  by  a  straight-line  interpolation 
between  calibration  points.  The  precision  in  reading  the  Hall  voltage 
was  about  ±0.02%.  The  value  of  the  magnetic  field  strength  was  known 
to  be  at  least  this  accurate,  and  the  accuracy  in  reading  the  voltage 
was  between  0.06  and  0.10%.  Accuracy  here  is  determined  by  the  ability 
to  take  repeated  measurements  of  the  same  value;  precision  is  deter¬ 
mined  by  the  number  of  places  to  which  this  measurement  may  be  read 
on  each  attempt. 


HALL  VOLTAGE  OUTPUT,  IN  MILLIVOLTS 


Rt-  SHUNT  RESISTOR 

P  -  LEEDS  a  NORTHRUP  POTENTIOMETER,  MODEL  K-3 
B  -  LARGE  CAPACITY  WET  CELL  BATTERY 
I  -  INPUT  CURRENT »  IOO  MA 


FIGURE  I.  Schematic  Diagram 


FIGURE  2.  Graph  of  Hall  Voltage  vs.  Magnetic  Field,  Showing 
Data  Points,  Best  Fit  Curve,  and  Best  Fit  Straight  Line 
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The  face  of  the  Hall  generator  was  estimated  to  be  within  3  degrees 
of  being  normal  to  the  field.  This  corresponds  to  an  error  of  0.1%  or 
less.  Error  due  to  temperature  change  and  instability  of  the  current 
was  negligible. 

Because  of  difficulties  in  fabrication,  the  leads  are  not  always 
attached  on  equipotential  lines.  Thus  one  may  read  a  voltage  across 
the  Hall  leads  when  the  input  current  is  flowing  even  though  the  field  is 
off.  This  may  also  effect  a  voltage  difference  when  the  field  is  non¬ 
zero  and  the  direction  of  the  field  is  reversed.  This  voltage  at  zero 
field  divided  by  the  input  current  is  defined  as  the  "resistive  zero  com¬ 
ponent."  For  the  device  reported  herein,  this  component  was  suffi¬ 
ciently  small  that  it  could  not  be  measured. 


SUMMARY  OF  RESULTS 


The  S  curve  deviated  from  a  linear  behavior  by  about  1%  with 
respect  to  the  value  at  10  kilogauss.  This  may  be  improved  by  a  factor 
of  3  by  using  the  sinusoidal  perturbation  mentioned  before.  Voltage 
readings  at  calibration  points  one  megacycle  apart  (adjacent)  differ  by 
less  than  0.15%;  hence  the  error  in  interpolating  between  adjacent  cali¬ 
bration  points  is  less  than  0.1%.  The  use  of  more  calibration  points 
and  better  accuracy  in  orienting  the  device  normal  to  the  field  would 
yield  an  error  estimated  at  less  than  0.05%.  Table  1  lists  the  data 
obtained  for  the  plot  shown  in  Figure  Z. 

Care  must  be  used,  especially  in  high  fields,  to  prevent  the  leads 
from  forming  large  loops  in  the  field.  Even  with  the  input  current  equal 
to  zero,  an  induced  voltage,  Uj,  may  be  set  up. 


U. 
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where  A  is  the  "inductive  zero  component"  expressed  in  cm^  and  dB/dt 
is  the  time  rate  of  change  of  the  field.  The  induced  current  may  easily 
be  high  enough  to  burn  the  leads  or  destroy  the  calibration  of  these  thin 
film  devices,  or  to  do  both. 


TABLE  1.  Calibration  Data  for  Hall  Generator* 


Probe 

Resonance  Frequency 
(±0.01%). 

Me 

Equivalent  Magnetic 
Field, 
kgauss 

Hall  Voltage, 

mv 

Hi 

10 

2.3487 

10.723 

11 

2.5835 

11.826 

12 

2.8184 

12.938 

13 

3.0533 

14.060 

14 

3.2881 

15.180 

15 

3.5280 

16.303 

16 

3. 7578 

17.430 

17 

3.9927 

18.563 

18 

4.2276 

19.696 

19 

4.4624 

20.852 

20 

4.6973 

21.998 

21 

4.9322 

23.154 

Li? 

14 

8.461 

40.366 

15 

9.065 

43.119 

16 

9.670 

46.003 

17 

10.274 

48.552 

* 

Type  FA-21,  Transverse  Field  Probe;  manufactured  by 
Siemens,  Munich,  Germany;  obtained  from  Instrument  Systems 
Corporation,  College  Point  56,  Long  Island,  New  York. 


SOLID  STATE  SPECTROSCOPY 


THE  INFRARED  ABSORPTION  SPECTRUM  OF  Ce3+  IN  LaCl3 

by 

R.  A.  Buchanan,  H.  H.  Caspers,  and  J.  Murphy 
INTRODUCTION 

In  a  previous  report*  the  problem  of  placing  a  trivalent  cerium 
ion  in  LaCl3  was  considered.  In  particular,  the  effect  of  the  crystal¬ 
line  field  symmetry  acting  on  the  free-ion  energy  levels  was  developed 
and  selection  rules  were  used  to  predict  the  number  of  absorption  lines 
due  to  the  electronic  transitions.  These  predictions  were  summarized 
in  Figures  2  and  3  of  the  report. 

The  spectrum  of  Ce^  +  in  LaCl3  in  the  spectral  region  from  1700  to 
4500  cm**  (2  to  6  microns)  has  now  been  obtained  and  will  be  reported. 
In  addition,  evaporated  films  of  LaCl3  have  been  prepared  and  their 
spectra  taken  in  the  infrared  region  between  60  and  700  cm"*.  A  quali¬ 
tative  discussion  of  these  spectra  is  included  here. 

The  infrared  spectrum  of  Ce^  +  in  CaF^  has  been  presented  by 
Feofilov”^  and  by  Kiss.^  Feofilov  has  noted  the  effects  of  charge  com¬ 
pensation  in  his  spectra,  and  Kiss  has  suggested  that  vibronic  lines 
are  present  in  his  spectrum.  No  attempt  was  made  to  analyze  these 
spectra  by  crystalline  field  theory,  although  Kiss  noted  a  crystalline 
field  splitting  of  about  1000  cm**.  Hutchison  and  Wong"*  reported  the 
paramagnetic  resonance  spectrum  of  Ce*'1'  in  LaCl3-  Their  work  indi¬ 
cates  which  state  is  the  ground  state.  This  is  essential  information  if 
predictions  of  the  absorption  spectrum  are  to  be  made. 


J.  Murphy,  H.  H.  Caspers,  and  R.  A.  Buchanan,  NAVWEPS 
Report  7237,  Quarterly  Report:  Foundational  Research  Projects,  July- 
September  1962,  pp.  111-134  (December  1962). 

2 

P.  P.  Feofilov,  Optics  and  Spectroscopy,  Vol.  6,  p.  150  (1959). 

3Z.  J.  Kiss,  Phys.  Rev.,  Vol.  127,  p.  718  (1962). 

4 

C.  A.  Hutchison,  Jr.  and  E.  Wong,  J.  Chem.  Phys.,  Vol.  29, 
p.  754  (1958).  “  ~ 
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EXPERIMENTAL  DETAILS 


Crystals  containing  0.1  and  5  mole-%  of  Ce*  +  in  LaCl3*  were  pre¬ 
pared  by  Professor  Eugene  Wong  of  the  UCLA  Physics  Department.  The 
preparation  technique  has  been  described  in  the  literature.'* 

The  crystals  were  ground  to  the  desired  shape  with  dry  emery  paper 
and  polished  on  a  cloth  lap  with  rouge  and  benzene.  The  polished  sam¬ 
ples  were  mounted  on  the  liquid  helium  Dewar  sample-holder  by  means 
of  spring  clips.  The  edges  of  the  sample  were  thickly  coated  with  sili¬ 
cone  vacuum  grease  to  provide  better  thermal  contact  between  the  crys¬ 
tal  and  sample-holder  at  the  lower  temperatures.  This  entire  operation 
was  carried  out  in  a  drybox. 

The  liquid  helium  Dewar  used  in  this  experiment  was  designed  spe¬ 
cifically  for  use  on  a  Cary  Model  14  spectrophotometer.  It  cools  the 
sample  by  conduction  of  heat  through  the  solid  copper  sample-holder. 
Since  the  sample  is  not  immersed  directly  in  the  liquid  helium,  the  true 
sample  temperature  is  unknown. 

The  spectra  were  recorded  in  the  near  infrared  region  by  using  the 
Cary  Model  14b  spectrophotometer  equipped  with  PbS  and  cooled  PbSe 
detectors  and  a  300  lines/mm  grating.  The  resolution  obtainable  with 
this  instrument  depends  upon  many  factors,  which  include  the  sample 
transmittance  and  the  wavelength.  In  the  unpolarized  spectra  being 
reported,  the  resolving  power  is  approximately  1000. 


DISCUSSION  OF  THE  SPECTRA 

Figure  1  shows  the  spectra  of  Ce^+  in  LaCl3  at  room  temperature, 
near  liquid  nitrogen  temperature,  and  near  liquid  helium  temperature. 
The  principal  effect  of  lowering  the  sample  temperature  is  to  narrow 
the  absorption  line  widths  and  thus  present  a  more  resolved  spectrum. 
The  room  temperature  spectrum  suggests  that  there  are  three  groups 
of  lines  centered  about  2250,  2900,  and  3400  cm"*.  Upon  being  cooled 
with  liquid  helium,  these  groups  of  lines  become  more  precisely  defined 
and  many  components  are  observed.  There  are  clearly  more  lines  here 
than  can  possibly  be  explained  by  electronic  transitions  alone.  The 


The  0.1  mole-%  sample  also  contains  about  0.1  mole-%  of  Gd-*+. 
Since  the  electronic  structure  of  Gd^  +  is  such  that  no  transitions  occur 
in  the  2-6  micron  region,  its  presence  will  be  ignored  in  the  following 
discussion.  It  should  be  noted  that  the  spectrum  of  the  crystal  that  con¬ 
tains  5  mole-%  Ce-*+  in  LaCl3  but  has  no  Gd^+  is  essentially  identical 
with  this  0.1  mole-%  Ce*+  sample. 


increase  in  absorption  in  the  1700  cm"*  (6  micron)  region  is  due  to  an 
instrumental  effect  and  is  not  characteristic  of  the  crystal  containing 
Ce3+  in  LaCl3. 

In  order  to  determine  the  electric  or  magnetic  dipole  nature  of  these 
transitions,  spectra  were  obtained  by  using  a  silver  chloride  plate  polar¬ 
izer.  These  spectra,  which  are  presented  in  Figure  2,  were  made  with 
the  sample  cooled  by  liquid  helium.  The  crystalline  c-axis  was  oriented 
perpendicular  to  the  slit  and  to  the  direction  of  light  propagation.  Effects 
of  polarization  are  present  in  these  spectra.  However,  an  unambiguous 
interpretation  of  these  effects  is  not  possible  at  this  time,  and  the  rea¬ 
sons  for  this  may  be  experimental.  The  spectrum  of  the  polarizer 
mounted  in  the  sample  compartment  in  such  a  way  that  the  S  was  per¬ 
pendicular  and  parallel  to  the  Blit  is  shown  in  Figure  A  large  polari¬ 
zation  of  the  spectrophotometer  optical  system  is  noted,  and,  in  addi¬ 
tion,  a  very  curious  crossing  of  the  curves  in  the  1-  to  2.5-micron  region 
raises  some  question  as  to  the  validity  of  any  polarization  spectra  taken 
on  this  machine.  Nevertheless,  the  spectra  of  the  polarizer  seem  well 
behaved  for  wavelengths  between  2.6  and  6.0  microns,  which  is  the  region 
of  primary  interest  in  this  investigation. 

The  crystalline  c-axis  was  located  by  using  crossed  polaroids.  The 
sample  was  mounted  so  that  its  c-axis  was  aligned  with  the  horizontal 
plane  and  was  perpendicular  to  the  direction  of  propagation.  The  polari¬ 
zer  was  then  mounted  externally  on  the  liquid  helium  Dewar  and  aligned 
to  be  parallel  and  perpendicular  to  the  horizontal  direction.  All  of  these 
alignments  were  made  by  eye  and  it  is  possible  that  some  misalignment 
occurred,  although  many  spectra  were  repeated  after  realignment  of 
the  polarizer  and  no  essential  changes  in  the  spectra  of  Figure  2  were 
observed. 

Perusal  of  the  lowest  temperature  spectrum  in  Figure  1  shows  that 
it  is  clearly  inconsistent  with  any  interpretation  based  only  on  electronic 
transitions  of  Ce^+  ions  on  a  single  symmetry  site.*  The  possibility  that 
Ce^+  ions  occupy  more  than  a  single  Ce  site  is  excluded  on  the  following 
grounds:  (1)  No  multiplicity  of  sites  has  been  reported  for  other  trivalent 
rare-earth  ions  in  LaCl3,  and  (2)  Hutchison  and  Wong^  have  studied  the 
ground  state  of  Ce^+  in  LaCl3  by  paramagnetic  resonance  techniques  and 
find  no  evidence  for  more  than  a  single  Ce^+  site. 

If  the  symmetry  of  the  La  site  were  not  C3h,  as  was  assumed  in 
the  previous  quarterly  report,*  the  number  of  multiplet  levels  in  the 
crystal,  and  hence  the  number  of  spectral  lines,  would  not  be  increased. 
The  minimum  degeneracy  of  the  multiplets  exists  for  C3h  symmetry 
since  each  level  consists  of  a  pair  of  Kramer's  conjugate  levels. 


76 


7 


FIGURE  3.  Transmission  of  a  AgCl  Plate  Polarizer  as  Measured 
with  the  Cary  Model  14b  Spectrophotometer 


If  only  sites  of  a  single  symmetry  are  available  to  the  Ce^+  ion  in 
the  LaClj  lattice,  it  would  seem  that  vibration  electronic  interactions 
are  the  only  plausible  explanation  for  the  "extra"  lines  observed  in  the 
spectrum.  Such  an  interaction  results  in  the  production  of  energy  levels 
that  are  separated  from  an  electronic  level  by  the  energy  of  a  lattice 
vibration.  Absorption  lines  associated  with  transitions  to  such  levels 
are  known  as  "vibronic"  lines.  The  problem  is  to  decide  correctly 
which  lines  of  the  spectrum  are  "pure  electronic"  and  which  ones  are 
"vibronic." 


The  spectra  of  a  crystal  of  LaClj  containing  5  mole-%  Ce^  +  are 
presented  in  Figure  4.  These  spectra  more  clearly  define  the  absorp¬ 
tion  lines  in  the  2250  cm-*  group.  Again,  the  principal  effect  of  lower¬ 
ing  the  sample  temperature  is  to  present  a  much  more  resolved  spec¬ 
trum.  In  addition  to  this,  the  intensities  of  the  absorption  lines  near 
2057  and  2100  cm-'  are  reduced  at  the  lower  temperature.  This  fact 
labels  these  absorptions  as  being  due  to  transitions  that  originate  on 
an  energy  level  above  the  ground  state. 


The  absorption  lines  at  2131  and  2168  cm“l  appear  to  be  narrower 
than  most  of  the  other  lines.  This  suggests  that  they  are  likely  candi¬ 
dates  for  assignment  as  pure  electronic  transitions.  In  support  of  this 
assignment,  it  is  to  be  noted  that  lines  are  observed  at  2210  and  2247 
cm-l.  These  lines  are  exactly  79  cm-*  above  the  assigned  electronic 
lines.  Such  lines  could  be  evidence  that  a  79  cm'*  lattice  vibration  is 
coupling  to  both  the  electronic  transitions.  Indeed,  if  one  searches 
for  37  separations  (2168  -  2131  =  37  cm'*)  in  the  region  between  2250 
and  2450  cm'*,  many  such  separations  can  be  found.  This  is  evidence 
that  many  other  lattice  vibrations  are  coupling  to  the  upper  electronic 
levels. 


The  possible  existence  of  vibronic  lines  in  the  spectrum  of  Ce^  + 
in  LaClj  has  stimulated  the  investigation  of  the  lattice  vibrations  of 
LaCl3.  A  complete  group  theoretical  analysis,  including  the  symmetry 
coordinates,  will  be  presented  in  the  following  paper  in  this  report. 

The  number  and  activity  of  the  normal  modes  of  fT  =  0  are  summarized 
in  Table  1.  Three  infrared  active  fundamental  absorption  bands  are 
predicted.  No  infrared  spectrum  of  LaCl3  has  previously  been  reported 
in  the  literature. 
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TABLE  1.  Selection  Rules  for  Transition 
From  the  Ground  State 


Operators 

Representation  of 
Modes 

Type  of  Spectrum 

(a  +  a  ),  a 

xx  yy  zz 

2  A  \ 

g  j 

a  ,  a 
yz  zx 

1  E.  > 
lg  ( 

Raman 

(a  -  a  ),  a 

xx  yy  xy 

3  E?  ) 

2g 

T 

z 

(T  ,  T  ) 
x  y 

2Elu$ 

Infrared 

2  B  \ 

g  / 

ZB  > 

U  | 

Inactive 

E2u 

LaCl3  FILMS 
Preparation 

LaC^,  particularly  in  thin  films,  reacts  rapidly  with  atmospheric 
H2O.  In  order  to  maintain  a  pure  LaCl3  sample,  it  was  found  necessary 
to  prepare  a  vacuum  cell  in  which  to  hold  the  evaporated  film.  The  cell 
mounted  in  the  evaporation  apparatus  is  shown  in  Figure  5. 

The  films  are  prepared  in  the  following  manner:  The  entire  vacuum 
cell  in  its  open  configuration  and  with  polyethylene  substrate  exposed  is 
placed  in  the  evaporation  chamber.  Next,  LaCl3  crystals  or  powders 
are  placed  in  the  platinum  boat  located  about  3  inches  directly  below  the 
substrate.  The  evaporation  chamber  is  evacuated  to  about  10^  torr 
and  the  boat  heated  sufficiently  to  cause  the  LaCl3  to  evaporate.  When 
a  film  of  the  desired  thickness  is  deposited  on  the  substrate,  the  current 
to  the  boat  is  turned  off  and  the  vacuum  cell  closed  and  sealed  by  means 
of  the  external  control  shaft,  which  enters  the  evaporation  chamber 
through  O-ring  seals.  Air  is  allowed  to  enter  the  evaporation  chamber, 
and  the  sealed  vacuum  cell  containing  the  evaporated  film  is  removed. 

The  spectrum  was  obtained  in  the  225  to  700  cm"'  region  on  a 
Beckman  IR-7  spectrophotometer  with  a  grating  and  Csl  foreprism 
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interchange  unit.  The  region  from  60  to  225  cm'*  was  covered  with  a 
Perkin- Elmer  98  monochromator  converted  for  grating  operation. 


The  spectra  of  two  evaporated  LaCl3  films  are  shown  in  Figure  6. 
These  are  preliminary  spectra  and  instrumentation  is  now  under  con¬ 
struction  that  will  permit  a  better  definition  and  hence  more  precise 
wave-number  assignment  of  the  absorption  bands. 

In  the  double-beam  IR-?  instrument,  the  vacuum  cell  with  the  sam¬ 
ple  was  placed  in  the  sample  beam  and  three  pieces  of  polyethylene  were 
placed  in  the  reference  beam  to  compensate  for  the  cell  windows  and 
substrate.  In  the  Perkin-Elmer  98  monochromator,  the  transmission 
signal  of  the  vacuum  cell  with  sample  was  divided  by  the  transmission 
signal  of  the  three  pieces  of  polyethylene  at  each  wavelength  to  deter¬ 
mine  the  percent  transmission  of  the  LaClj  film  at  that  particular  wave¬ 
length.  The  crosses  on  Figure  6  represent  the  percent  transmission  of 
the  three  pieces  of  0.040-inch  polyethylene. 

The  three  absorption  bands  near  136,  165,  and  212  cm'*  are  appar¬ 
ently  associated  with  the  three  infrared  active  £T  =  0  modes  (see  Table  1). 
The  evaporated  films  consist  of  randomly  oriented  crystallites  and,  there¬ 
fore,  polarization  data  are  not  available.  Without  polarization  data,  any 
assignments  of  the  Eju  and  Au  modes  to  the  observed  absorption  bands 
are  uncertain.  Assignments  of  the  Eju  and  Au  modes  may  be  possible 
by  comparing  the  far  infrared  spectra  of  the  isomorphous  compounds 
.LaBrj  and  UCI3  with  that  of  LaCl3.  This  work  is  now  in  progress. 

It  is  of  considerable  interest  to  examine  the  vibronic  spectrum  of 
Ce^+  to  see  if  these  infrared  active  lattice  frequencies  couple  to  the 
electronic  levels.  Table  2  shows  that  the  136  and  212  cm'*  modes  are 
present  in  the  vibronic  spectrum.  The  165  cm'*  mode  does  not  yield 
such  good  agreement.  If  a  mode  of  173  cm'*  is  assumed,  good  agree¬ 
ment  is  obtained. 


TABLE  2.  Vibronic  Spectrum  of  Ce^ 


Electronic 
cm"  1 

Vibration 
+  cm'*  = 

Predicted  Vibronic 
cm"* 

Observed  Vibronic 
cm'  * 

2131 

136 

2267 

2268 

2131 

173 

2304 

2302 

2131 

212 

2343 

2342 

2168 

136 

2304 

2302 

2168 

173 

2341 

2342 

2168 

212 

2380 

2381 
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Several  explanations  for  the  failure  of  the  165  cm"'  mode  to  be 
observed  in  the  vibronic  spectrum  can  be  presented.  Until  better  far- 
infrared  data  actually  confirm  or  reject  the  wave-number  value  165 
cm'*,  speculation  on  this  matter  will  be  withheld. 

CONCLUSIONS 

From  these  observations,  it  must  be  concluded  that  there  is  strong 
evidence  for  the  assumption  that  lattice  frequencies  of  the  LaCl3  host 
are  coupling  to  electronic  levels  of  the  Ce^  +  ions.  As  a  result  of  this 
coupling,  many  more  lines  appear  in  the  absorption  spectrum  in  the  4- 
to  5-micron  region  than  would  otherwise  exist  if  only  "pure"  electronic 
transitions  were  involved.  No  attempt  has  been  made  to  interpret  the 
groups  of  lines  near  2900  and  3400  cm~^.  These  lines  are  too  far  re¬ 
moved  from  the  2131  and  2168  cm*^  lines  to  be  vibronic  lines  associated 
with  these  electronic  transitions.  They  may  be  vibronic  lines  of  other 
electronic  transitions,  however. 

It  is  tempting  to  associate  some  of  the  lines  near  2900  and  3400 
cm"*  with  vibration-induced  electronic  transitions.  The  selection  rules 
of  Murphy  et  al.^  were  based  on  the  static  crystalline  field  model  and 
can  be  shown  to  break  down  in  the  vibrating  crystalline  field.  This  as¬ 
signment  leads  to  a  larger  crystalline  field  splitting  than  is  generally 
observed  for  rare-earth  ions  in  LaClj.  It  is,  however,  consistent  with 
the  1000  cm'^  crystalline  field  splitting  observed  by  Kiss^  for  Ce^  in 
CaF2. 

The  possibility  that  impurities  are  present  in  the  spectra  has  been 
considered  but  it  is  felt  that  this  possibility  cannot  be  the  explanation  of 
the  additional  lines,  for  the  following  reasons:  (1)  The  spectral  regions 
in  the  shorter  wavelength  infrared  and  visible  have  been  examined  and 
contain  no  absorption  bands.  This  would  seem  to  eliminate  another 
rare-earth  ion  as  the  impurity,  (2)  The  spectrum  of  Ce^+  in  LaF3 
shows  a  spectrum  similar  to  Ce3+  in  LaCl3.  This  would  seem  to  imply 
that  the  spectra  presented  are  characteristic  of  Ce^+.  The  possibility 
of  some  other  impurity,  such  as  an  organic  molecule,  cannot  be  posi¬ 
tively  eliminated  at  this  time. 

The  Ce^+  spectrum  seems  to  defy  interpretation  in  all  host  lattices 
in  which  it  has  been  studied.  It  is  necessary  to  obtain  the  positions  of 
the  electronic  levels  and  establish  the  crystalline  field  parameters 
before  any  quantitative  understanding  of  the  Ce^+  spectrum  can  be 
obtained. 
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SYMMETRY  COORDINATES  FOR  LaCl3 
by 

J.  Murphy  and  H.  H.  Caspers 


INTRODUCTION 

For  processes  in  which  only  a  single  phonon  is  involved,  the  maxi¬ 
mum  intensity  of  an  absorption  line  in  the  vibrational  spectrum  of  a 
lattice  corresponds  to  a  fiT  =  0  mode.  This  comes  about  because  a  IT  =  0 
infrared  active  mode  has  associated  with  it  a  large  induced  electric  di¬ 
pole  moment  in  comparison  with  a  IT  /  0  mode  of  the  same  branch.  For 
a  similar  reason,  the  peaks  in  the  Raman  spectrum  are  associated  with 
IT  =  0  vibrations.  There  seems  to  be  good  agreement  between  the  absorp¬ 
tion  peaks  found  in  the  infrared  lattice  spectrum  and  the  "vibronic"  lines 
observed  in  the  optical  absorption  spectra  of  rare-earth  ions  in 

crystals. 

The  factoring  of  the  secular  determinant  for  the  frequencies  of  the 
IT  =  0  modes  is  greatly  facilitated  by  the  introduction  of  symmetry  co¬ 
ordinates.  These  symmetry  coordinates  reduce  the  determinant  to 
block  diagonal  form.  The  subdeterminants  for  the  frequencies  now 
involve  only  those  symmetry  coordinates  that  transform  as  the  same 
row  of  the  same  irreducible  representation  of  the  factor  group.  It  is 
desirable  to  establish  a  procedure  for  determining  sets  of  symmetry 
coordinates  that  form  bases  in  terms  of  which  the  normal  modes  of  the 
lattice  can  be  expressed.  The  linear  combination  of  symmetry  coordi¬ 
nates  which  form  a  normal  mode  depend  upon  the  selection  of  significant 
force  constants.  Naturally,  an  exact  solution  would  require  an  infinite 
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number  of  force  constants  for  an  infinite  lattice.  However,  it  is  often  a 
good  approximation  to  consider  the  interactions  only  up  to  a  few  nearest 
neighbors.  One  tries  to  express  all  of  the  vibrational  frequencies  in 
terms  of  as  small  a  number  of  force  constants  as  possible.  As  one  will 
discover  in  applying  this  formalism,  certain  lattices  prove  to  be  rather 
difficult  because  of  either  too  many  particles  in  the  primitive  cell  or  too 
low  a  symmetry,  or  both.  There  are,  however,  a  large  number  of  lat¬ 
tices  for  which  the  symmetry  coordinates  can  readily  be  reduced. 


LATTICE  VIBRATIONS  IN  LaCl3 

LaCl3  crystallizes  in  a  hexagonal  dipyramidal  type  of  symmetry 
and  belongs  to  the  space  group  C63/m  (C(yy?’)J  Figure  7  shows  the  par¬ 
ticles  in  the  primitive  cell.  (The  x,y,z  coordinate  system  indicated  will 
be  used  in  labeling  particle  displacements  in  Table  9.)  According  to 
Zachariasen,?  the  dimensions  of  the  primitive  cell  are  a3  =  a£  =  7.468  A, 
a3  =  4.366  A.  The  particle  positions  are  expressed  as  fractions  of  the 
dimensions  of  the  primitive  cell  and  are 

2  La:  (1/3,  2/3,  1/4),  (2/3,  1/3,  3/4) 

6 Cl:  (a,  b,  1/4),  (1  -  b,  a  -  b,  1/4),  (1  -  a  +  b,  1  -  a.  l/4) 

(a  -  b,  a,  3/4),  (1  -  a,  1  -  b,  3/4),  (b.  1  -  a  +  b,  3/4) 

where  a  =  0.375  and  b  =  0.292.  The  origin  of  coordinates  is  chosen  as 
the  point  marked  0.  Figure  8,  which  is  taken  from  International  Tables 
for  X-Ray  Crystallography, 8  shows  all  of  the  symmetries  of  the  primi¬ 
tive  cell.  The  symmetry  operations  that  leave  Figure  8  invariant  are 
listed  in  Table  3.  The  k  =  0  factor  group  of  Cfch^  can  be  seen  to  be  iso¬ 
morphic  with  the  underlying  point  group  C^;  therefore,  these  two  groups 
have  the  same  character  table.  From  Figures  7  and  8,  one  can  deter¬ 
mine  that  the  point  symmetry  at  the  site  of  a  La  ion  is  C3j1,  while  the 
point  symmetry  at  a  chlorine  site  is  Cs,  i.  e.,  (E,  tr^).  Table  4  lists 
the  irreducible  characters  of  C(,h^>  which  is  isomorphic  with  the  point 
Cfch-  Tables  5  and  6  list,  respectively,  the  irreducible  characters  of 
C3h  and  Cs.  The  groups  C3h  and  Cs  are  subgroups  of  C^h- 

It  is  of  interest  to  determine  the  transformation  properties  of  these 
modes  at  the  sites  of  the  La  and  Cl  ions.  A  mode  with  given  transforma¬ 
tion  properties  can  only  occur  if  the  particles  that  make  up  the  lattice 
can  move  and  maintain  those  transformation  properties.  The  point  group 

- 
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FIGURE  7.  Positions  of  Ions  in  the  Primitive  Cell. 
The  large  circles  represent  lanthanum  ions  and  the 
smaller  circles,  chlorine  ions.  Ions  in  the  1/4  plane 
are  shown  filled  in;  those  in  the  3/4  plane  are  open. 
Zero  (0)  is  the  origin  of  coordinates  for  the  ion  posi¬ 
tions  given  in  the  text.  The  particle  displacements 
listed  later  in  Table  9  are  written  with  respect  to  the 
indicated  x,y,z  coordinate  system  where  the  z  axis  is 
out  of  the  plane  of  the  paper. 
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LEGEND: 


Q  -60*  rotation  followed  by  inversion 

£  -60*  rotation  followed  by  translation 

of  aj/2  where  is  the  z  dimen¬ 
sion  of  the  primitive  cell 

0  -180*  rotation  followed  by  a$/2  trans¬ 

lation 


^T,T  the 


planes  at  l/4  a3  and  3/4  a3  are 
eflection  planes 


inversion 


FIGURE  8.  Symmetry  Operations  in  the  Cfch^  Primitive 
Cell.  All  rotations  are  about  axes  passing  through  the 
symbol  and  perpendicular  to  the  plane  of  the  paper. 
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TABLE  3. 


Symbol 

U|Tn) 

{C3|0},  {C32  1 0 } 

{S3 |0) 

{ S35  1 0 >  =  {l|0}{C6|0} 


h 1 0  > 
{I|0} 


{s6|0},  {s65|0} 
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2 

Operations  of  the  Space  Group  C^h 


Description 

Lattice  translation:  Tn  =  nj  Iti  +  n£a2  +  n3^3 
where  | a  1  |  =  |  a2  |  t  I&33  | 

Three-fold  rotation  about  axes  parallel  to  the 
c-axis  through  symbols  and  £ 

Three -fold  rotatory  reflection  about  axis  paral¬ 
lel  to  crystalline  c-axis  passing  through  @ 

Horizontal  reflection  about  plane  perpendicular 
to  c-axis  passing  through  l/4  and  3/4  planes 

Inversion  through  centers  indicated  by  circles  o 

Screw  diad  normal  to  paper  passing  through 
symbols  £ 


Screw  hexads  normal  to  paper  passing  through 
symbols  £ 

Six-fold  rotatory  reflections  with  axes  parallel 
to  c-axis  and  passing  through  symbols  ^ 


TABLE  5.  Character  Table  for  C3h 


C3h 

E  C3 

S' 

=35 

% 

S3 

A1 

1  1 

1 

1 

1 

1 

tr  • 

n  -e* 

-e 

-£ 

1 

-e  * 

U  -e 

-€  * 

-£* 

1 

-£ 

A" 

l  l 

1 

-1 

-1 

-1 

|PII 

n  -e  * 

-£ 

e 

-1 

€* 

ll  -£ 

-e  * 

€  * 

-1 

£ 

Note:  e  =  exp  (2rri/6) 

TABLE  6. 

Character  Table  for  C8 

C 

E 

cr . 

s 

h 

A' 

1 

1 

A" 

1 

-1 

at  the  site  of  each  particle  must  be  a  subgroup  of  the  lattice  space  group. 
Each  representation  of  the  point  group  with  which  the  factor  group  is 
iBomorphic  (let  us  call  it  G)  can  be  expressed  as  linear  combinations  of 
the  representations  of  the  point  group  at  each  particle  site.  In  other 
words,  the  representations  of  G  become  reducible  representations  of 
the  point  groups  at  the  particle  sites.  Let  us  suppose  that  the  particle 
coordinates  (or  displacements)  transform  according  to  the  representa¬ 
tions  yx,  Yy>  and  Yz  °f  the  point  group  gp,  where  gp  denotes  the  point 
group  at  the  particle  site.  Upon  reducing  the  representation  Tj  of  G  to 
those  of  gp,  at  least  one  of  the  Yx*  Yy>  or  Yz  must  appear  at  least  once 
in  the  reduction;  otherwise  no  mode  of  the  lattice  that  transforms  as  Tj 
will  exist  in  which  the  particle  at  the  gp  symmetry  site  will  move.  On 
the  other  hand,  suppose  at  the  La  site  the  representation  T j  reduces  to, 
say,  Ya  +  Yb  +  Yz*  where  a  and  b  are  not  x  or  y.  In  the  lattice  mode 
that  transforms  as  rj,  the  La  ion  can  only  move  in  the  z  direction  since 
only  the  z  coordinate  has  the  right  transformation  properties  for  that 
mode.  This  situation  requires,  then,  that  if  a  representation  of  the 
space  group  does  not  contain  any  representation  according  to  which 
the  particle  coordinates  transform  for  any  type  of  particles  of  the  crys¬ 
tal,  there  is  no  mode  of  the  lattice  that  transforms  as  Using  these 

arguments,  one  can  draw  a  correlation  chart1^  (Figure  9),  which  shows 

o 
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FIGURE  9.  Correlation  Chart  Showing  the  Reduction  of  C^  to  C3h 
and  C3  Symmetry.  The  symbol  (Tx,  Ty)  means  that  x  and  y  trans¬ 
lations  transform  into  linear  combinations  of  one  another  according 
to  the  representation  next  to  which  they  appear.  The  symbols  (Tx), 
(Ty),  and  (Tz)  mean  that  the  x,  y,  and  z  translations  transform  into 
themselves  times  a  constant. 


the  reduction  of  the  space  group  representations  into  the  representations 
of  the  particle  site  symmetries,  and  from  it  determine  which  particles 
move  in  the  various  modes. 

One  can,  in  fact,  determine  the  number  of  unit  cell  modes  and  their 
transformation  properties  quite  apart  from  the  usual  method  of  reducing 
the  reducible  representation  of  the  primitive  cell  into  the  irreducible 
representations  of  the  space  group.  These  two  methods  serve  as  a  use¬ 
ful  check  on  each  other.  From  Figure  9,  we  note  first  of  all  that  modes 
of  vibration  (including  translations  of  the  whole  crystal)  are  allowed  for 
every  representation  of  C(,h-  Secondly,  in  modes  that  transform  as  Ag, 
Bu,  Ejg,  and  E2u*  on^Y  the  chlorine  ions  can  take  part;  the  lanthanum 
ions  do  not  move. 

Given  a  set  of  n  independent  coordinates,  one  can  write  n  linearly 
independent  combinations  of  these  coordinates.  This  implies  that  there 
are  n  linearly  independent  modes  of  vibrations  including  complete  trans¬ 
lations  and  rotations.  Therefore,  for  each  degree  of  freedom  we  can 
say  there  is  one  mode.  In  a  symmetry  coordinate,  the  particle  displace¬ 
ments  that  transform  as  an  even  representation  are  such  that  two  par¬ 
ticles  which  are  related  by  inversion  through  an  inversion  point  in  the 
lattice  must  be  moving  in  opposite  directions,  whereas  the  displace¬ 
ments  of  two  such  particles  in  an  odd  mode  are  in  the  same  direction. 

If  there  is  a  particle  sitting  at  an  inversion  point  of  the  primitive  cell, 
an  even  mode  cannot  involve  the  displacements  of  that  particle  since  it 
cannot  be  going  in  two  directions  at  once.  In  such  a  case,  the  number 
of  even  and  odd  modes  would  not  be  equal.  In  LaC^,  however,  there 
are  no  particles  at  the  inversion  points.  In  Figure  9,  an  odd  representa¬ 
tion  is  listed  under  the  even  representation  to  which  it  is  related  by  a 
change  in  sign  under  inversion,  i.  e.,  if  in  an  E2g  symmetry  coordinate 
the  particle  displacements  on  one  side  of  the  inversion  center  are 
changed  in  phase  by  180  deg,  the  resulting  set  of  displacements  would 
transform  as  an  Eju  representation.  Hence,  if  we  can  determine  all  of 
the  even  symmetry  coordinates,  we  can  find  all  of  the  odd  ones,  and 
vice  versa. 

An  Ag  mode  (Figure  9)  involves  only  the  x,y  motion  of  the  chlorines. 
In  LaClj  each  chlorine  ion  can  be  transformed  into  any  other  chlorine 
ion  by  some  operation  of  the  space  group.  Under  such  a  transformation, 
the  x  and  y  components  of  the  displaced  ion  become  some  linear  combi¬ 
nation  of  the  x  and  y  coordinates  of  the  ion  into  whose  position  it  has 
been  carried  by  the  symmetry  operation.  This  two-dimensional  trans¬ 
formation  implies  two,  and  only  two,  independent  solutions  or  eigen¬ 
functions  of  the  transformation  matrices.  Hence  there  are  two  Ag 
modes.  The  number  of  Ag  modes  and  Bu  modes  must  be  equal  since 
the  Bu  modes  can  be  obtained  from  the  Ag  modes  by  changing  the  sign 
of  the  particle  displacements  in  the  lower  half  of  Figure  7.  We  have 


thus  used  up  four  linear  combinations  of  the  chlorine  x,y  coordinates, 
or,  as  far  as  numbers  are  concerned,  four  out  of  the  twelve  coordinates. 
The  E2g  and  Eju  modes  (which  are  equal  in  number)  involve  the  displace¬ 
ments  of  both  the  chlorine  and  the  lanthanum  ions  in  the  x,y  plane.  There 
are  two  lanthanum  ions  in  the  primitive  cell  and,  therefore,  four  new 
x,y  coordinates.  This  means  that  we  once  again  have  twelve  coordinates, 
and  so  twelve  linear  combinations  or  modes  must  be  divided  equally 
among  the  &Zg  arK*  Equ  representations.  There  are,  therefore,  3  E^g 
modes  and  3  Eiu  modes  [(3  +  3)  x  2  =  IE]  .  For  the  representations  at 
the  lower  half  of  Figure  9,  the  two  linear  combinations  of  lanthanum  z 
displacements  must  be  divided  equally  between  the  Au  and  the  Bg  repre¬ 
sentations.  The  six  chlorine  z  displacements  are  divided  as  follows: 


—  •  *  .  E 

Eu 

A  .  . 
u 

.  1 

.  .  ■  •  2 

B  .  . 

.  1 

g 

Each  E  representation  must  have  two  linearly  independent  modes  since 
an  E  representation  is  two-dimensional.  The  remaining  two  modes  are 
divided  between  the  Au  and  the  Bg  representations.  These  results,  which 
are  tabulated  in  Table  7,  are  the  same  as  those  that  would  be  determined 
by  the  usual  group  theoretical  method  of  reducing  the  reducible  character 
of  the  unit  cell. 9 


TABLE  7.  Numbers  and  Species  of  Modes  in  LaClj 


Represen¬ 

tation 

j 

No.  of  Represen¬ 
tations 

No.  of 

Modes 

Types  of 
Ions 
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Polari¬ 
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A 

g 

E 
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Cl 

(x,y) 

B 

u 

E 

2 

Cl 

(x.y) 

E2g 

3 

6 

La  +  Cl 

(x.y) 

Elu 

3 

6 

La  +  Cl 

(x,y) 

A 
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E 

2 

La  +  Cl 

Z 

B 

g 

E 

2 

La  +  Cl 

z 

Elg 

1 

E 

Cl 

z 

EEu 

1 
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Cl 
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Such  a  reducible  character  iB  shown  in  Table  8  and  is  obtained  by 
performing  each  symmetry  operation  of  the  space  group  on  the  particles 
in  the  primitive  cell  and  finding  how  many  such  particles  transform  into 
themselves  or  into  particles  that  are  primitive  translations  away.  The 
set  of  these  numbers  multiplied  by  the  traces  of  the  corresponding  trans¬ 
formation  matrices  is  the  character  of  the  reducible  representation  and 
is  equal  to  the  sum  of  the  irreducible  representations  of  Cf>h  weighted 
by  the  numbers  shown  in  the  second  column  of  Table  7. 


TABLE  8.  Character  Table  for  the  Reducible  Representation 
of  the  Primitive  Cell  Modes 


X(R)  24  00  0  0  00  -4  0  80  -4 


From  the  correlation  chart,  we  have  determined  not  only  the  num¬ 
bers  and  types  of  normal  modes,  but  also  which  particles  move  and 
whether  the  motion  is  in  the  x,y  plane  or  in  the  z  direction,  This  infor¬ 
mation  allows  us  to  find  a  set  of  symmetry  coordinates — basis  functions 
from  which  normal  modes  are  determined.  For  example,  we  have  deter¬ 
mined  that  there  are  three  pairs  of  modes  that  transform  according  to 
the  E2g  representation  in  LaClj.  In  each  of  these  modes,  both  the  La 
and  Cl  ions  are  moving.  Furthermore,  the  motion  is  restricted  to  be 
parallel  to  the  x,y  plane.  Three  modes  transform  as  the  first  row  of 
the  representation  E^g  and  three  transform  as  the  second  row.  Any  set 
of  displacements  of  the  particles  that  transform  as  the  first  row  of  E2g 
must  be  expressible  as  a  linear  combination  of  normal  modes  which 
transform  as  that  first  row.  Conversely,  a  normal  mode  is  expressible 
as  a  linear  combination  of  linearly  independent  sets  of  displacements 
that  transform  according  to  the  representation  of  interest.  Conse¬ 
quently,  if  we  can  find  three  linearly  independent  sets  of  displacements 
that  transform  as  the  first  row  of  the  representation  E2g.  then  the  three 
normal  modes  will  be  expressible  as  linear  combinations  of  these  three 
"symmetry  coordinates." 

The  normal  modes  are  the  linear  combinations  of  symmetry  coordi¬ 
nates  that  completely  diagonalize  the  secular  determinant.  We  will  not 
go  so  far  as  to  calculate  the  secular  determinant  in  this  paper  as  we  are 
only  interested  in  demonstrating  the  method  of  obtaining  symmetry 
coordinates. 
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DETERMINATION  OF  SYMMETRY  COORDINATES 

A  symmetry  coordinate  S,  which  transforms  according  to  a  one¬ 
dimensional  irreducible  representation  of  the  factor  group,  has  the  fol¬ 
lowing  property: 


RS  =  S’  =  kS  (1) 

where  R  represents  an  element  of  the  group  and  k  is  a  constant  factor  of 
unit  modules.  In  the  case  of  an  abelian  group  such  as  the  factor  k 

can  be  complex.  In  a  non-abelian  group,  k  is  restricted  to  ±1.  In  the 
case  of  complex  k's,  the  representations  that  are  complex  conjugates  of 
one  another,  i.e.,  factors  k  and  k*,  are  generally  taken  together  and 
the  combination  is  treated  as  a  two-dimensional  representation  of  the 
group.  In  the  present  study,  it  is  convenient  to  maintain  distinct  one¬ 
dimensional  representation  of  the  C(,h  point  group. 

As  an  example  of  the  procedure  to  be  followed,  let  us  choose  the 
representation  listed  in  Table  4  as  E-Zg  >  that  is,  the  first  row  of  E^g. 

Consider  the  operation  {C& |  a 3 / 2 }  of  Table  3  about  the  point  C  of 
Figure  7.  This  operation  has  the  following  effect  on  the  particles: 

1- 5  7  -  2'  =  2 

2— 8  8  -  3 1  =3 

3— 6  6  —  4'  =  4 

4- 7  5  -  1 '  =  1 

The  primes  indicate  particles  that  are  removed  from  the  particles  with¬ 
in  the  unit  cell  by  a  primitive  translation.  Since  we  are  considering 
k  =  0  modes,  the  primed  and  unprimed  particles  move  exactly  in  phase. 
This  has  been  indicated  by  writing  2'  =  2,  etc.  According  to  equation  (1) 
and  the  character  table  for  the  representation  E^g*,  the  new  sets  of  dis¬ 
placements  are  related  to  the  old  by 


In  other  words,  the  transformed  particle  2  is  leading  the  untransformed 
particle  8  by  120  deg.  A  second  application  of  {C^|a3/2}  or,  equiva¬ 
lently,  C  3  |  0 } ,  carries  particle  2  into  particle  3  (or  an  equivalent  posi¬ 
tion  displaced  from  the  first  by  a  primitive  translation).  We  see  from 
the  Cfcii  character  table  and  Figure  7  that  the  transformed  particle  2  is 
lagging  the  untransformed  particle  3  by  120  deg.  According  to  the  defi¬ 
nition  given  here  of  a  symmetry  coordinate,  we  are  at  liberty  to  choose 
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the  displacement  of  any  particle  arbitrarily  and  then  operate  on  this 
particle  with  the  group  elements  to  find  the  displacements  of  the  other 
particles.  It  is  convenient  to  choose  the  motion  of  our  starting  particle 
(say,  particle  2  in  Figure  7)  to  be  circular.  By  carrying  out  the  pro¬ 
cedure  outlined  above,  we  arrive  at  Figures  10a,  10b,  and  10c.  Except 
for  modes  that  involve  translations  of  the  whole  crystal,  the  lanthanum 
motions  relative  to  the  chlorine  motions  are  chosen  so  that  the  linear 
momentum  of  the  primitive  cell  is  zero.  It  can  easily  be  verified  that 
there  are  three,  and  only  three,  symmetry  coordinates  (barring  multi¬ 
plication  of  each  particle  motion  by  a  constant  phase  factor)  that  involve 
circular  motion  of  all  of  the  particles  and  transform  as  EZg  •  The  impli¬ 
cation  is  not  that  this  is  the  only  set  of  symmetry  coordinates  which 
transform  this  way.  One  can  form  an  infinite  number  of  sets  of  three 
linearly  independent  motions  of  the  three  sets  shown.  The  symmetry 
coordinates  that  transform  as  E2g^  are  the  complex  conjugates  of  those 
that  transform  as  E2g^  and  are  shown  in  Figures  11a,  lib,  and  11c.  It 
is  not  convenient  for  purposes  of  calculating  the  vibrational  frequencies 
to  use  the  symmetry  coordinates  just  derived.  Linear  combinations  that 
are  the  sum  and  difference  of  the  complex  conjugate  pairs  turn  out  to  be 
better,  since  in  these  all  particles  pass  through  their  equilibrium  posi¬ 
tions  simultaneously.  Furthermore,  the  new  motions  are  rectilinear. 

The  sum  and  difference  of  Figures  10a  and  11a  are  shown  in  Figures  12a 
and  13a  respectively.  Figures  12b  and  13b  are  the  sum  and  difference 
of  Figures  10b  and  lib,  while  Figures  12c  and  13c  are  the  sum  and  dif¬ 
ference  of  Figures  10c  and  11c. 

Except  that  there  is  only  one  coordinate  for  each  particle,  the  analy¬ 
sis  of  the  z-polarized  symmetry  coordinates  is  carried  out  in  the  same 
fashion  as  the  x,y  motions.  The  displacement  of  one  of  the  particles  is 
chosen  arbitrarily  and  the  motions  of  all  of  the  particles  are  determined 
with  respect  to  it.  In  the  z-polarized  symmetry  coordinates  that  trans¬ 
form  as  Ejg  and  E2U,  the  three  chlorine  ions  in  each  triangle  must  be 
120  deg  out  of  phase  with  one  another.  Once  again,  however,  the  sums 
and  differences  of  the  symmetry  coordinates  so  derived  are  more  con¬ 
venient  for  calculational  purposes. 

The  symmetry  coordinate  for  the  kT  =  0  modes  of  LaCl3  are  listed 
in  Table  9.  Those  shown  that  transform  as  E  representations  are  actu¬ 
ally  the  sums  and  differences  of  the  symmetry  coordinates  first  arrived 
at,  as  discussed  above. 

The  particle  displacements  are  listed  as  (i,j,k)  where  i,  j,  and  k 
are  the  components  of  the  displacement  vectors  located  on  the  particles 
and  written  with  respect  to  the  x,y,z  coordinate  system  shown  in  Fig¬ 
ure  7.  Some  of  the  lanthanum  displacements  are  multiplied  by  a  factor 
which  is  required  in  order  to  reduce  the  linear  momentum  of  the  unit 
cell  to  zero. 


FIGURE  11.  Symmetry  Coordinates  That  Transform  as 


FIGURE  12.  Symmetry  Coordinates  That  Transform  as  + 

(a,  b,  and  c  appear  in  Table  9  as  S5a,  Ssy,,  and  S5e,  respectively.  ) 


FIGURE  13.  Symmetry  Coordinates  That  Transform  as  E^gl  *  ^Zg% 
(a,  b,  and  c  appear  in  Table  9  as  Sfca,  S6b.  and  S6c.  respectively.  ) 
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where  ro  =  chlorine  mui  end  M  *  lanthanum  mtn. 


SPECTRAL  EMITTANCE  OF  SOLIDS 


SPECTRAL  EMITTANCE  MEASUREMENTS  OF 
SOME  COMMERCIAL  TRANSPARENT  SOLIDS 

by 

D.  L.  Stierwalt,  D.  D.  Kirk,  and  J.  B.  Bernstein 

A  continuing  study!  is  being  conducted  at  this  Laboratory  on  the 
spectral  emittance  of  solids.  The  theory,  techniques,  and  experimental 
apparatus  used  in  acquiring  the  spectral  emittance  of  opaque  and  trans¬ 
parent  solids  in  the  50-200*C  temperature  range  in  the  infrared  region 
of  the  spectrum  have  been  discussed  in  detail  in  a  previous  report. ^  In 
this  report  are  presented  data  on  two  commercial  transparent  solids 
that  have  some  applications  in  the  near  infrared  region. 

Two  samples  of  transparent  solids,  1  mm  thick,  were  made  avail¬ 
able  through  the  courtesy  of  the  Eastman  Kodak  Company.  Termed 
"Irtran  3"  and  "Irtran  4,"  they  consist  of,  respectively,  polycrystalline 
calcium  fluoride  and  polycrystalline  zinc  selenide.  Figure  1  shows  the 
spectral  emittance  of  Irtran  3  at  100*C  and  Figure  2  shows  the  emittance 
of  Irtran  4. 


^NAVWEPS  Report  7237,  Quarterly  Report:  Foundational  Research 
Projects,  July-September  1962,  p.  135  (December  1962). 

2 

NAVWEPS  Report  5981  (NOLC  487),  Quarterly  Report;  Foundational 
Research  Projects,  Qctober-December  1959,  p.  41  (April  I960). 
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R.  L.  Conger,  Code  42  1 

C.  A.  Roberts,  Code  422  1 

K.  I.  Crowder,  Code  43  1 

R.  F.  Potter,  Code  43  1 

D.  L.  Stierwalt,  Code  43  1 

J.  B.  Bernstein,  Code  432  1 

D.  D.  Kirk,  Code  432  1 

R.  A.  Buchanan,  Code  433  1 

H.  H.  Caspers,  Code  433  1 

J.  Murphy,  Code  433  1 

H.  Piller,  Code  433  1 

C.  P.  Haber,  Code  44  1 

W.  S.  Harris,  Code  441  1 

R.  E.  Panzer,  Code  441  1 

W.  C.  Spindler,  Code  441  1 

C.  M.  Douglas,  Code  442  1 

D.  L.  Herring,  Code  442  1 

F.  C.  Essig,  Code  45  1 

J.  R.  Alday,  Code  452  1 

E.  W.  Seeley,  Code  452  1 

R.  W.  Yancey,  Code  452  1 

B.  F.  Husten,  Code  50  1 

L.  L.  Parker,  Code  50  1 

H.  A.  Bulgerin,  Code  51  1 

R.  L.  Bauer,  Code  53  1 

J.  J.  Nastronero,  Code  54  1 

R.  J.  Hardy,  Code  57  1 

F.  C.  Alpers,  Code  71  1 

W.  F.  Meggers,  Jr.,  Code  72  1 

T.  B.  Jackson,  Code  73  1 

A.  C.  Sevy,  Code  75  1 

Technical  Library,  Code  234  2 
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